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ABSTRACT 

We present a physics-based statistical theory of a force-free magnetic field in the corona above a turbulent 
accretion disk. The field is represented by a statistical ensemble of loops tied to the disk. Each loop evolves 
under several physical processes: Keplerian shear, turbulent random walk of the disk footpoints, and recon- 
nection with other loops. To build a statistical description, we introduce the distribution function of loops over 
their sizes and construct a kinetic equation that governs its evolution. This loop kinetic equation is formally 
analogous to Boltzmann's kinetic equation, with loop-loop reconnection described by a binary collision inte- 
gral. A dimensionless parameter is introduced to scale the (unknown) overall rate of reconnection relative to 
Keplerian shear. After solving for the loop distribution function numerically, we calculate self-consistently 
the distribution of the mean magnetic pressure and dissipation rate with height, and the equilibrium shapes of 
loops of different sizes. We also compute the energy and torque associated with a given loop, as well as the total 
magnetic energy and torque in the corona. We explore the dependence of these quantities on the reconnection 
parameter and find that they can be greatly enhanced if reconnection between loops is suppressed. 
Subject headings: 



1. INTRODUCTION 

Power-law components in the X-ray spectra of accreting 
black holes are attributed to hot, tenuous, quasi-spherical 
plasmas called accretion disk coronae (hereafter ADCe) 
by analogy with the solar corona, though comptoniza- 
tion rather than atomic lines or bremsstrahlung is thought 
to be the main emission mechanism in accreting system s 
dBisnovatvi-Kogan & Blinnikovlll976t iLiang & Pried fT977h . 
The vertical extent of ADCe is open to doubt since they are 
spatially unresolved and since the electron temperature in- 
ferred from high-energy spectral cutoffs, T e ~ lOOkeV, is 
typically small compared to the virial temperature of the 
ions, Tj ~ lOOMeV. Direct simulations of magnetorotational 
(hereafter MRI) turbulence in radiation-pressure-dominated 
disks show very la rge density contrasts and bulk velocities 
dTurner et al.1 120021) . so perhaps the po wer-laws are made 
within the disk itself (ISocrates et al.ll2004l) . 

If a comptonizing corona indeed exists, then it is difficult 
to avoid the conclusion that it should be magnetically dom- 
inated. The electrons th emselves cannot store much energy 
dMerloni & Fabian 2001): their Compton cooling time is at 
most comparable to the local dynamical time at luminosities 
> 10~ 2 LEdd and radii < 20GMbh/c 2 . Ions at virial tempera- 
tures would store much more energy but could no t transfer it 
efficie ntly to the electrons by Coulomb collisions dRees et al.1 
1982). Also, as in the solar case, magnetic fields are probably 
needed to convert mechanical energy of the optically thick 
regions into coronal heat. And since accretion is believed 
to be driven by MRI turbulence within the disk proper, it is 
expected that fields float up into the corona by Parker and 
interc hange instabilities (Gal eev et al.1 Il979t iTout & P ringle 
fl992h . 

Despite widespread recognition of these points, modelers 
of ADC emission rarely concern themselves with the dy- 
namics of coronal magnetic fields. This is perceived to be 
too hard; certainly the long and arduous struggle to under- 
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stand the heating of solar corona — based on much more abun- 
dant data — tends to discourage hopes of solving the corre- 
sponding problem for accreting system s any time soon (e.g. 
IWalsh & Irelandll2003l: lKlimchukl l2006). Attempts to model 
coronae through direct th ree-dimensional (3D) magnetohy- 
drodynamic shearing-box (Miller & Stone 2000; Hiros e et all 
2006) and global dMachida et al.l2000l) simulations have been 
made, however. To date, such simulations suggest that while a 
magnetically dominated region does form, its vertical extent, 
when defined by distributions of shear stress or dissipation 
rate, exceeds that of the optically thick regions only modestly. 

There are nevertheless good reasons to expect dynami- 
cally dominant coronae and to question the contrary evidence 
from simulations. On the one hand, a strongly magnetized 
corona may be required to e xtract power from black hole spin 
(Blandford & Znaiekj [T977|) or from plunging gas ins ide the 
marginally stable orbit (Gammi 3 [1999 ; Krolik 1999b . or to 
drive a wind from the disk (Blandford & Payne 1982). A less 
widespread motivation, which however we feel strongly, is to 
reduce self-gravity in the outer parts of accretion disks in ac- 
tive galactic nuclei (at r > IQr'GM^/c 1 ): a magnetized corona 
or wind might transport angular momentum more quickly 
than stresses limited to the optically thick layer, and thereby 
redu ce the mass dens ity within the disk for a given accretion 
rate (IGoodmanll2003l) . 

On the other hand, codes designed for pressure-dominated 
plasmas may not be reliable when applied to magnetically 
dominated coronae. Shearing-box calculations cannot be 
trusted to represent magnetic structures much larger than a 
disk scale height (H) unless all three dimensions of the box 
are ^> H, which has not yet been achieved. Global simula- 
tions are unable to resolve thin disks and extended coronae 
unless the grid is made coarser in the corona than near the 
midplane, which increases numerical diffusion in the corona. 
This is not serious when H ~ r, as for the innermost parts 
of near-Eddington or radiatively inefficient accretion disks, 
where the X-ray evidence for coronae is strongest; but it 
would be a severe limitation for simulations of the marginally 
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self-gravitating parts of AGN disks, where H/r < 10~ 2 . These 
problems of spatial dynamic range will eventually be over- 
come with sufficient computer power. 

A more fundamental difficulty for simulations has to do 
with magnetic reconnection. We will argue that the efficiency 
of reconnection 1 is crucial to the storage and dissipation of 
magneti c energy in ADCe, as it appears to be in the solar 
corona dParkerll 19831 119881: Hjjmchuk 2006). Reconnection 
in solar flares is observed to be "fast," meaning that annihilat- 
ing field lines come together at a significant fraction of the lo- 
cal Alfven speed. This is not well understood; because of the 
very high conductivity of the corona, MHD predic ts reconnec- 
tion rates slo wer by many orders of magnitude (ISweetlll958l: 
lParker|[T957l) . Collisionless plasma effe cts outside conve n- 
tional resistive MHD may be necessary (Uzdensky 2 007alfbl 
and references therein). Another mystery is what triggers fast 
reconnection, which does not always occur immediately but 
seems to require special c onditions that have not been fully 
identified (Klimchuk 2006). The overall rate of magnetic dis- 
sipation depends both on the rate at which reconnection events 
are triggered and on the speed of reconnection during such 
events. 

These complexities are elided by the astrophysical MHD 
codes used to study magnetorotational turbulence. Often in 
these codes, reconnection is purely numerical, that is, pro- 
duced by truncation errors due to limited grid resolution. Ex- 
plicit diffusivities are sometimes used, but if these are large 
enough make truncation errors unimportant, then they are 
necessarily orders of magnitude larger than astrophysical re- 
ality. Until recently, the general view among MRI simula- 
tors seems to have been that the microphysics of reconnec- 
tion is unimportant. This view may be inspired by analogies 
with the dissipation of kinetic energy in hydrodynamic tur- 
bulence, both compressible and incompressible. Supersonic 
turbulence involves shocks, and as is well known, shock dis- 
sipation is independent of transport coefficients in the limit 
that these are small. On the other hand, three-dimensional in- 
compressible turbulence involves inertial cascades such that 
dissipation, though occurring on small (viscous) scales, is en- 
tirely controlled by the dynamics on large scales and proceeds 
at rates that are again asymptotically independent of transport 
coefficients. 

Whether or when turbulent magnetic dissipation is sim- 
ilarly independent of small — and therefore numerically 
unresolved — scales is an open question. Recent work in- 
dicates that magnetorotational turbulenc e is sensitive to the 
magnetic Prandtl number P m = vjr\ dLesur & Lon garetti 
2007; Fromang et al. 2007), at least in the range of fluid and 
magnetic Reynolds numbers (Re = VL/v, Re m = P m Re) ac- 
cessible to direct simulations. It may be that Re„, and Re be- 
come unimportant when both are sufficiently large, but this 
has not yet been established even for kinematic dynamos, 
where the fie ld is dynamically unimportant on all scales by 
const ruction dBoldyrev & Cattan"eo1l2004t ISchekochihin et alJ 
2007, and references therein), much less for MRI turbulence. 
Even if this is true of systems in which fluid motions dominate 
the energy density, as they are presumed to do near the mid- 
plane of an accretion disk, the answer could be different for 
magnetically dominated systems such as ADCe. Perhaps rel- 

1 It would be natural to write "rate of reconnection" here, but that this 
phrase is often used to mean the speed at which field lines of opposing polar- 
ities approach one another in a localized reconnection event, whereas we are 
concerned here with the global efficiency of reconnection in reducing mag- 
netic energy. 



evant here is the case of turbulence in the presence of a dom- 
inant mean field, which has becom e somewhat better under- 
stood since the seminal paper of Goldr eich & Sridharl (1 19951 
hereafter GS95). Cascades do exist in such turbulence, with 
different scalings along and perpendicular to the mean field, 
so that the large scales are insulated from details at the re- 
sistive and viscous scales. However, there are at least two 
important differences between such systems and the nearly 
force-free ADCe contemplated in this paper. In the former 
systems, the Alfvenic propagation time along the field is as- 
sumed to be longer than the timescale of the cascade, so that 
the turbulent dynamics are local; by contrast, communication 
along the entire length of a line-tied coronal loop is effectively 
instantaneous compared to the timescale on which energy is 
injected into the loop by footpoint motions. The magnetic 
dynamics are therefore nonlocal, and not appropriately char- 
acterized as turbulent; a better description is a progression of 
force-fre e equilibria driven by ch anges at the boundary. [See, 
however, Rapp azzo et al.1 (120071) for a contrary view.] Sec- 
ondly, in GS95's Alfvenic cascades, the deviations from the 
mean field are small, and the cascade does not affect the en- 
ergy in the large-scale mean field. In this paper, by contrast, 
we are concerned with coronal flux loops that reconnect with 
one another at large angles between their field lines. A sig- 
nificant fraction of the loops' magnetic energy may be liber- 
ated in such reconnection events, or in the relaxation to a new 
force-free equilibrium following topological changes brought 
about by reconnection. 

It might be hoped that the MRI simulations could predict 
the total dissipation rate of the corona, if not its vertical dis- 
tribution, because the rate of injection of energy to the corona 
is determined at its base, where thermal and kinetic energies 
dominate. This is a false hope, however. The rate of work 
done on the corona by the disk is proportional to the mag- 
netic stress tensor at its base, specifically the rz and <f>z com- 
ponents of the stress. Insofar as the corona is approximately 
force-free, its total energy is expressible in terms of a bound- 
ary integral involving the same stress components. Thus if 
the coronal energy and magnetic configuration are incorrectly 
calculated, then the coronal dissipation rate is probably also 
incorrect. More concretely, because most of the kinetic en- 
ergy available from the disk is the large-scale differential ro- 
tation rather than local turbulence, coronal field loops with 
large radial separations between their footpoints may be par- 
ticularly important for the energy input; such loops are not 
possible in shearing-box simulations whose radial dimensions 
are no larger than the disk scale height, and even if the dimen- 
sions were increased, spurious reconnection might suppress 
the large loops. In short, for geometrical reasons and because 
they do not model reconnection correctly, present-day MRI 
simulations may underestimate the energies and dissipation 
rates of disk coronae. 

The dynamics of ADCe may not be fully understood with- 
out great progress on all of the fronts described above: more 
powerful computations, better understanding of fast recon- 
nection, and of course more incisive observations. Since all 
of this may take years or decades to accomplish, our purpose 
here is to try to imagine, in a disciplined way, some aspects of 
th at ultimate understandin g. Our approach is clearly indebted 
to iTout & Pringlel d 19961 henceforth TP96) but is richer in 
physical elements. We model the ADC as an assembly of 
closed magnetic loops with footpoints on the disk; open field 
is probably important but is deferred to a later paper because 
we do not wish to address winds here. We do not resolve 
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the dynamics within the optically thick disk at all but treat 
the base of the corona as a dynamic boundary. It is assumed 
that the disk thickness is much less than its radius in the parts 
of the ADC that we model and that the loop lengths (L) lie 
at intermediate scales r). Thus the lower bound- 

ary is conceived as an infinite plane. Small new loops are 
injected from this boundary, and existing loops are energized 
by the Keplerian shear if their footpoints lie at different radii. 
Loops reconnect in pairs according to prescribed rules with a 
frequency scaled relative to the shear by a dimensionless pa- 
rameter. Each reconnection results in a new pair involving 
the same four footpoints differently connected, rather than a 
single loop as in TP96. The loop population is described by a 
distribution function over the length and orientation of the dis- 
placement from the negative to the positive-polarity footpoint 
(rather than the length alone as in TP96), and the processes 
of injection, stretching, and reconnection are incorporated in 
an integro-differential kinetic equation for the evolution of the 
loop distribution function. 

Our model is not appropria te for all forms of ADCe. Fol- 
lowing iLiang & Pried dl977l) . we presuppose a "sandwich" 
geometry in which the optically thick but geometrically thin 
disk coexists with its corona at the same radius. Such a ge- 
ometry seems most appropriate to high/soft and very high 
states of galactic X-ray sources where the X-ray continuum 
shows both thermal and power-law features, and also to cases, 
both galactic and extragalactic, where a relativistically broad 
iron Ka line indicates that an optically thick disk, fluoresc- 
ing under illumination by hard X-rays from the corona, ex- 
tends down to the marginally stable orbit or beyond (e.g. 
IWilms et al] 1200 U IMiller et all |2002|) . In low/hard states 
where only a power-law is seen, it may be that the inner parts 
of the thin disk have evap orated so that those regions are all 
"corona" (lEsin et alJll997l) : alternatively, the thin disk may 
persist down to the marginally stable orbit, but the corona may 
take the form of a mildly relativistic outflowin g wind or jet, 
whose emission is direct ed away from the disk (iBeloborodovl 
119991: IMiller et"a l. 2006). The coronal model presented here 
would not apply to either of these cases without substantial 
modification. 

In this paper, we explore what controls the magnetic energy, 
integrated stress, and dissipation rate of the corona; specifi- 
cally, we explore the roles of shear and reconnection in this 
balance. Our overarching motivation is to determine under 
what conditions the corona contributes importantly to outward 
angular momentum transport: that is, to the torque that drives 
accretion through the disk. Secondary goals are to examine 
the distribution of energy, stress, and dissipation with height. 
We recognize that because of a number of questionable as- 
sumptions and simplifications, our model will hardly be the 
last word on this subject. We hope , however, at least to set up 
a target for future simulations to aim at, and to draw attention 
to certain quantities that could be extracted from existing sim- 
ulations, such as the rate of emergence of flux dipoles (small 
loops) from the disk. 

The outline of our paper is as follows. § [2] introduces our 
conception of magnetic loops and the loop distribution func- 
tion. §[3]explores some properties of loops in equilibrium with 
the mean-field pressure exerted by neighboring loops: shape, 
maximum height, and energy. §|4]constructs the kinetic equa- 
tion, including the rules for reconnection and other processes 
important to the evolution of the loop distribution. The nu- 
merical set-up used to solve the loop kinetic equation, includ- 
ing the boundary conditions at both small and large scales, is 



described in §[5j In the same section we also present numeri- 
cal solutions and discuss their implications for ADCe. In §|6] 
we discuss limitations of our model in the light of these re- 
sults, and we indicate ways in which the model might be made 
more realistic. We attempt to relate what we have done to the 
present theoretical understanding of ADCe, and we discuss 
how simulations in the near future might be used to calibrate 
some features of our model, for example the rate of emer- 
gence of small loops from the disk. Finally, § [7] summarizes 
our main conclusions. 

2. STATISTICAL DESCRIPTION OF THE ADC MAGNETIC FIELD 

As noted in §1, there are fundamental physical similarities 
between the formation processes of the solar and accretion 
disk coronae. At the same time, one has to keep in mind the 
important differences between them. These include: (^dif- 
ferences in the underlying sub-photospheric turbulence (ther- 
mal convection vs. MRI); (2) a strong large-scale differential 
rotation and shear in disks, whereas small scales dominate 
the shear in the Sun; (3) the possibility of strong magneto- 
centrifugally- and radiation-driven winds from an accretion 
disk, as compared to the relatively weak thermally driven so- 
lar wind; (4) Compton cooling of the disk corona in black hole 
systems; (5) a greater separation of spatial scales in thin disks 
between the disk thickness and its radius, as opposed to the 
solar case, where the convection zone spans a significant frac- 
tion of the Sun's radius, which results in the generation and 
emergence of very large magnetic structures associated with 
sunspots in active regions. Finally, we have only one Sun, 
whereas there is a large variety of astrophysical accretion-disk 
systems, some of which have coronae. 

In addition to the above, there is an important basic dif- 
ference in our observational capabilities: whereas the Sun is 
so close that we can spatially resolve individual events and 
structures in the solar corona (such as flares, loops, etc), such 
resolution is not available for ADCe. Therefore, we can study 
only the spatially-integrated spectral and timing properties of 
the disk corona. This fact provides a strong motivation for 
focusing on a statistical description for the magnetic field in 
the ADC. 

Here is, briefly, the basic physical picture of the magne- 
tized corona above a turbulent accretion disk. The corona 
is a dynamic, self-organized system th at can be represented 
by a statistical ense mble of flux loops dTout & PringldfT996l : 
Hughe s"et alJl2003h . The loops continuously emerge out of 
(and submerge into) the disk as a result of magnetic buoyancy. 
Once above the surface, they constantly evolve due to a num- 
ber of physical processes. They are twisted and stretched by 
the differential keplerian rotation and by the random motions 
of their footpoints on the disk's surface, which causes the in- 
dividual loops to inflate. As a result, the magnetic field in the 
corona becomes non-potential and highly stressed; an appre- 
ciable amount of free magnetic energy can thus be stored in 
the corona. However, in the process of twisting and expan- 
sion the loops may undergo internal disruptions due to MHD 
instabilities and also may reconnect with other loops. Such 
relaxation events manifest themselves as flares; they bring 
the field closer to the potential state and thus enable the in- 
flation process to resume. At the same time, reconnection 
between loops sometimes produces more spatially-extended 
magnetic structures in the corona (the coronal "inverse cas- 
cade"). Finally, all these complicated processes occur repeat- 
edly over and over, simultaneously on various spatial scales. 
Thus, the corona can be viewed as a boiling magnetic foam, 



FIG. 1 . — Schematic view of ADC as an ensemble of many magnetic loops. 



FIG. 2. — Closed magnetic loop as the main structural element of a magne- 
tized accretion disk corona. 



in which magnetic loops repeatedly swell and grow because 
of the magnetic energy pumped into them by the footpoint 
motions and then snap and contract back due to reconnective 
disruptions and sometimes merge to form bigger structures. 

What is the appropriate statistical language for describing a 
chaotic, highly intermittent magnetized corona above a turbu- 
lent accretion disk? At a most general level, we describe the 
corona by an ensemble of some elementary magnetic struc- 
tures. Each of these fundamental individual constituents of 
the corona is characterized by a small set of primary phys- 
ical parameters. The individual elements evolve with time 
according to certain rules that reflect the relevant physical 
processes that we believe are most important in shaping the 
corona. These processes generally include various forms of 
interaction between elements. Mathematically, the evolution 
rules are represented by stochastic (Langevin) equations of 
motion of magnetic elements in the primary parameter space. 
Since we are interested in a statistical description, we intro- 
duce a distribution function of our magnetic elements in the 
primary-parameter phase-space. Correspondingly, one of our 
main goals is to derive the equation for the evolution of this 
distribution function, using the equations of motion of indi- 
vidual elements. This is done by analogy with the way the 
Boltzmann kinetic equation for the particle distribution func- 
tion in a gas is derived in Statistical Mechanics, but obviously 
is more ad hoc in our case. Finally, there are several impor- 
tant integral quantities in our theoretical framework, which 
are related to moments of the distribution function. These 
self-consistent quantities represent the mean-field interaction 
between the magnetic elements and they affect the evolution 
of the distribution function. 

2.1. The Loop Distribution Function 

Based on the above discussion, we shall now build a calcu- 
lable model of the corona. Our first task is to select the most 
appropriate and most fundamental elementary magnetic con- 
stituents of the corona. We shall then need to select the most 
natural set of parameters describing these elements. 

Guided by the analogy with the solar corona, we shall 
use simple (anchored in the disk at both ends, see Fig. [2]i 
magnetic loops, or flux tubes, as our fundamental magnetic 
elements — the main structural constituents of the corona. 
This choice is influenc ed by the existing t heoretical work in 
both solar physics (e.g . IHughes et al.ll2003l) . and also in ADC 
(Tout & Pringle 1996). Such loops represent the closed mag- 
netic field corresponding to zero net vertical flux through the 
disk. This is a natural assumption for the case when the mag- 
netic field in the corona comes from the flux emergence of the 



field generated by the dynamo in the disk itself. In this paper 
we shall assume that this is indeed the case. In principle, how- 
ever, one may also wish to consider a more general situation 
where, in addition to the closed coronal loops, there is also 
an external large-scale open magnetic field through the disk, 
such as may be coming from the central star or the interstellar 
medium (see § 16. 11 1. 

For simplicity, we shall characterize each loop by only 
two primary parameters: (i) the radial footpoint separation 
Ar = r_ — r+, and ( ii) the azimuthal footpoint separation r A<f> = 
Ay = y- —y+ (see Fig. O. Thus, we shall measure Ar and Ay 
from the "+" magnetic footpoint to the "-" magnetic foot- 
point. This means that Ar and Ay can be positive or nega- 
tive depending on the orientation of the loop. Alternatively, 
sometimes we will use an equivalent representation in terms 
of the loop's projected length (the distance between the loop's 
footpoints): L= (Ar 2 + Ay 2 ) 1 / 2 ; and the orientation angle, 
6, measured clockwise with respect to the toroidal direction: 
6 = arctan(Ar/Ay). 

In a more general description, one may enlarge the parame- 
ter space to include additional parameters, such as the mag- 
netic flux A^ contained within a loop, or the loop's twist 
(see § 16.21 ). However, since our goal here is to build the sim- 
plest version of this already very complicated theory, we shall 
assume that all the loops have the same magnetic flux A\l/ 
and, furthermore, that they are not twisted. The latter assump- 
tion means that the magnetic field within each loop is purely 
potential, i.e., that the bulk of the corona is nearly current-free 
and all the coronal currents flow along inter-loop boundaries. 

In addition to the above two primary parameters, we shall 
also need some secondary parameters describing a given loop, 
such as the loop's overall shape, its maximum height Z top 
above the disk; and its thickness at a given height, d(z). These 
quantities will be useful for estimating the loop expansion rate 
and for analyzing binary interaction (reconnection) of loops 
with each other. In our model, these secondary parameters 
are uniquely determined by the primary ones in combination 
with the self-consistent mean field B(z) (see § 13.21 ). 

Note that a loop carrying a finite flux A^ has a certain fi- 
nite thickness (in the radial and azimuthal directions) at the 
disk surface. For a typical loop, this thickness is generally of 
order H. Therefore, instead of a pair of footpoints that a field 
line would have, a finite-thickness loop has a pair of footspots. 
Thus, we need to be a little bit more precise in our definitions 
of Ar and Ay. We shall define them as the radial and az- 
imuthal separations between the centers of the two footspots, 
that is between the two footpoints of the central field line of 
the loop (magnetic axis for a twisted loop). 
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Following iTout & Pringld (119961) . we introduce the 
loop distribution function, F(Ar, Ay), defined so that 
F( Ar, Ay)dArdAy is the number of loops with the values of 
primary parameters in the range ([Ar, Ar + dAr], [Ay, Ay + 
dAy]), per unit disk area. Alternatively, we may write the dis- 
tribution function in terms of loop length L = (Ar 2 + Ay 2 ) 1 / 2 
and orientation 9 = arctan(Ar/Ay), i.e., F(L,8). The overall 
normalization of the distribution function is determined by the 
requirement that loops cover the entire disk surface; it will be 
discussed in more detail in § 13.21 

2.2. Role of Magnetic Reconnection in ADCe 

At the most basic level, the corona (either solar or ADC) 
is perfectly conducting almost everywhere. However, as it 
evolves driven by the complex turbulent motions of the mag- 
netic footpoints on the surface, the corona may develop nu- 
merou s current sheets on a variety of scales ( Parker! [1972L 
1983). These current sheets are possible sites of dissipation 
of magnetic energy via reconnection. In fact, reconnection 
is one of the most essential nontrivial physical processes that 
govern the complex dynamical behavior of the corona. In par- 
ticular, it controls the vertical extent of the corona {e.g., the 
coronal magnetic scale-height Hg). Indeed, if reconnection 
were too efficient, then the coronal field would be nearly po- 
tential and Hb ~ H; then, the free magnetic energy stored in 
the corona would be small, as would the magnetic dissipation 
rate. 

On the other hand, if no reconnection were allowed at 
all, then, magnetic loops would, over time, grow bigger and 
bigger in height because of the differential Keplerian rota- 
tion. Unable to dissipate, magnetic energy would continu- 
ously accumulate in the corona as the characteristic magnetic 
scale-height Hg increases, essentially linearly in time. This 
phase would continue until Hg becomes comparable to the 
disk radius, Hg ~ r ^> H. After that, radial gradients be- 
come important and the subsequent evolution would enter a 
qualitatively different regime characterized by accelerated ex- 
pansion of the magnetic loops, which effec tively would be- 
come open, perhaps even in a finite time Cy an Balle gooiieri 
19911: iLvnd en-Bell & Boily] [1991: IaTvI 11995b ISturrocketaT. 
19951 IUzdenskyll2002l) . The corona would then consist of a 
dense forest of open flux tubes of alternating polarity sepa- 
rated by a multitude of current sheets. Unless there are sig- 
nificant mass-loaded winds (violating the force-free assump- 
tion), the power pumped into the corona would then go down, 
and the accumulated free magnetic energy would saturate at 
a value cor responding to a fully open (split-monopole) mag- 
netic field (I Alyll 1 99 it ISturrockli 1991h . Although this asymp- 
totic energy would be very large, of order r/H larger than that 
of the fully closed potential field, it would still remain finite; 
this is because the toroidal magnetic field at the disk surface 
and hence the work done on the coronal magnetic field by 
the Keplerian disk shear would both go to zero. Similarly, 
the angular momentum exchange between different parts of 
the disk due to coronal loops would also go down. As new 
flux tubes emerged from the disk, the magnetic forest would 
become ever more dense. The corona would thus look very 
different from what we expect. We thus see that reconnection 
is necessary for maintaining a meaningful statistical steady 
state. It enables open field lines to close back and thus restores 
the magnetic connection between different parts of the disk. 
This, in turn, facilitates angular-momen tum transport vi a the 
coro nal magnetic field [coronal MRI. (lGoodmanll2003l) : see 
also dHevvaerts & Priesi)ll989HPavlidou et alj|2001l) 1. Strong 



magnetic dissipation and large torque thus require some in- 
termediate reconnection efficiency, neither so rapid as to keep 
the field nearly potential, nor so slow as to allow it to become 
fully open; in both limits, the torque and energy dissipation 
rate vanish. 

Another reason why reconnection is important is that a 
growing magnetic loop may reconnect with another one con- 
nected to a very different place on the disk. This pro- 
cess may lead to an "inverse cascade " of magnetic loops 
(Tout & Pr ingld [T996I) . It is an important avenue towards 
building up a population of loops with large radial footpoint 
separation. Indeed, whereas Keplerian differential rotation in- 
creases the azimuthal footpoint separation of a loop, it does 
not affect its radial footpoint separation. Therefore, without 
reconnection, the radial footpoint separation of a coronal flux 
loop would change only relatively slowly by the random walk 
of its footpoints due to the underlying disk turbulence. In 
principle, the footpoints will drift radially apart in direct re- 
sponse to the angular momentum transfer by the coronal loop 
itself (we call this process "coronal MRI"). The characteris- 
tic velocity of this drift is on the order of B 2 /47rS57, where 
£ is the surface density of the disk. The resulting relative in- 
crease in Ar on the rotation-period timescale is of the order 
of 5Ar/Ar ~ fi 2 /47r£Arf> 2 . Using £ ~ Hp and H ~ c s /Q, 
where c s and p are the sound speed and the gas density within 
the disk, we can estimate that SAr/Ar ~ (Va/c s ) 2 H / Ar = 
f3~ l H/Ar, where V% = B 2 /4irp is the Alfven speed within 
the disk. Thus, since we are mostly interested in large loops, 
Ar ^> H, we see that Ar cannot grow appreciably without re- 
connection. Thus reconnection is necessary for the coronal 
"inverse cascade ". 

In addition, magnetic reconnection in the corona regulates 
the fraction of the magnetic flux that is open at any given time 
and also the e ffective radial transport of a large-scale v ertical 
magnetic field ISpruit & Uzdenskyl (120051) : iFiskl (120051) . Both 
of these processes are important for establishing large-scale 
disk outflows. 

Finally, as in the solar corona, reconnection is believed to be 
the main mechanism of releasing the accumulated magnetic 
energy, leading to coronal heating and observed high-energy 
coronal emission. 

3. THE SELF-CONSISTENT CORONA 

In principle, the loop distribution function F(L, 9) contains 
enough information to fully describe the statistical magnetic 
structure of the corona. This means that, once F(L,9) is 
known, one should be able to answer most of the questions 
posed in §Q] In particular, one should be able to derive the 
actual shapes and heights, Z top (L), of coronal loops, the dis- 
tribution of magnetic energy with height, B 2 (z)/8ir, the en- 
ergy £(L) associated with a loop of a given size, the torque 
transmitted by the coronal magnetic field, etc. In this section 
we demonstrate how to do all this. 

3.1. Equilibrium Shape of a Loop in a Stratified Atmosphere 

First, we shall work out the correct shape of an isolated 
slender (with a cross-sectional diameter d <C L) untwisted 
loop A carrying magnetic flux A "J, immersed in a medium 
with some isotropic but, in general, nonuniform pressure P(z) 
(see Fig. [3). This pressure represents the confining magnetic 
pressure of all other loops; thus, for actual calculations, it will 
be convenient to write P(z) as P(z) = B 2 {z)/%n. The shape of 
the loop is then determined by the requirement that the loop 
be in magnetostatic equilibrium with this external pressure. 
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FIG. 3. — Untwisted isolated loop confined by external isotropic pres- 
sure P(z). 



First, the local pressure balance across the loop gives us the 
magnetic field strength inside it as a function of height: 



B(z)=^UP(z) = B(z). 



(1) 



Then, since the magnetic flux is constant along the loop, 
we can write the local cross-sectional area a(z) of the loop in 
terms of B(z)'- 

A* 

a(z) = wy (2) 

Now let us discuss equilibrium shape x(z) of a slender loop 
as a whole, represented by the shape of the loop's central field 
line. First, we would like to note that, for a curved but un- 
twisted loop confined by a uniform external pressure Pq, it is 
impossible to find an equilibrium shape. Indeed, since the 
field inside the loop has no twist and is purely axial (i.e., runs 
along the loop), it is a potential field produced by perpendicu- 
lar (to the direction of the loop) currents flowing on the loop's 
surface. Therefore, at any given location along the loop, the 
magnetic field strength is slightly nonuniform in the cross- 
loop direction: it drops off as 1 /R, where R is the local curva- 
ture radius ("major radius" in tokamak terminology). In other 
words, the magnetic force balance between the magnetic pres- 
sure and the magnetic tension inside a curved loop means that 
there must be a magnetic pressure gradient to balance the ten- 
sion force due to the curvature. Therefore, the magnetic field 
on the underside (the "inboard", in tokamak terminology) of 
the loop is larger than that on the upside (the "outboard"). On 
the other hand, however, the magnetic field at each point on 
the surface of the loop has to be in pressure balance with the 
external pressure P(z). If P(z) = Po = const, this pressure bal- 
ance means that the magnetic field has to be uniform along 
the loop's boundary. Thus, we have a clear incompatibility 
between the assumption that the loop is curved but untwisted 
and the condition of local force balance with an external uni- 
form pressure. 

This phenomenon can also be understood using the notion 
of Pfirsch-Schluter currents, a well-known concept in toka- 
mak physics. Since the magnetic field inside an untwisted 
flux tube is strictly in the axial direction (along the tube), the 
currents that produce it flow on the skin of the tube in the 
perpendicular direction. But since the currents have to close, 
and since the tube is curved, the surface current density (i.e., 
the current per unit length along the tube) is bigger on the 
underside than on the upside. Then, according to Ampere's 
law, the magnetic field is also stronger on the underside, and 
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FIG. 4. — Infinitesimal loop segment. 



thus cannot be in pressure equilibrium with a uniform ex- 
ternal pressure on both sides simultaneously. The only way 
a curved loop can be in equilibrium with uniform external 
pressure, is when the surface current density is also uniform. 
Since the current has to be conserved, this requires that some 
surface current should also flow along the loop (the so-called 
Pfirsch-Schluter current), and hence the magnetic field must 
be twisted. 

These considerations show that an untwisted flux tube con- 
fined by a uniform external pressure has to be straight. If, 
however, the pressure is not uniform, then an equilibrium 
shape for a curved tube can be found, as we now show. 

Let us consider a case with pressure P(z) decreasing mono- 
tonically with height. Consider a slender loop in the (x,z) 
plane, symmetric with respect to x = 0, with its two footpoints 
at x = ±L/2 (see Fig. [3]). The shape of the loop as a whole, 
described by a symmetric function z(x), is determined by the 
perpendicular (to the magnetic field) force-balance between 
tension and the pressure-gradient forces. The force-balance is 
established at every point along the loop. To analyze it at a 
given location on the loop, let us consider a small loop seg- 
ment. We shall denote the arc-length along the loop, measured 
from its left footpoint, by /. For clarity of discussion, let us 
represent the loop segment locally by a slightly curved cylin- 
der of length 51 <§C L, cross-sectional area a(l), and curvature 
radius R(l) (see Fig. Let us also introduce the angle a(l) 
between the magnetic field and the vertical direction: 



dz 

cos a = — . 
dl 



(3) 



For definiteness, we consider the ascending leg of the loop, 
where B z > 0, so that < a < tt/2. The local curvature ra- 
dius R(l) of the loop is related to a(l) via 



1 

R 



da 

~dl 



(4) 



Here we treat B(z) as a known function and our goal is to de- 
termine the geometrical shape of the loop described by a(z). 

The magnetic tension force on the loop segment acts in the 
direction of its curvature radius R and is equal to 

B 2 

Vtemion = -r^a(0^, (5) 

The projection of the external-pressure force on the loop 
segment onto R can be written as 

dP 



bfvp = - 



dz 



sin a a(l)Sl . 



(6) 
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Then, the force-balance condition can be written simply as 



— = - sin a 



dP 

dz' 



dB 



B 2 1 
4^ R 

or, making use of equation (|T), 

B(z) 

R dz 

Combining this equation with the geometrical relation (0), 
we immediately obtain: 

da I d log B 
~dl~~~~ 



- sin a - 



R 



dz 



(7) 



(8) 



(9) 



[Note that, as one can immediately see, a magnetic loop can 
be in equilibrium with a uniform [B(z) = const] external pres- 
sure only if it is straight, a = const.] 
Using ([3J, we get 



cot a 



da dlogB 



dz 



dz 



sin a(z) = 



C 

Wzj 



(10) 



It is interesting to note that the same result can be obtained 
in a simple and elegant way by using a variational principle, 
namely, by by minimizing the loop's magnetic energy E magD = 
(A*/87r) J B(l)dl (see §133) viewed as a functional of z(x). 
Using the relationship dl = dx sin -1 a, we have 



SE m , sn [z(x)] = ^-5 J B[z(x)] Vl+z' 2 (x)dx 



A* 
~8T 



B'(z)-B(z) 



l+zf- 



Sz(x)dx 



.(id 



From the condition 5E mlign = 0, we thus immediately get 
B = CVl+z' 2 (x) = C dl C 



dx sin a 



(12) 



which, taking into account that B = B(z), is the same as the 
above result ( fTOb . 

The result ( fTOb means that the horizontal (jc) component of 
the magnetic field is constant along the loop: 



B x (l) = B sin a = B(z) sin a = C = const . 



(13) 



The integration constant C is just equal to the magnetic 
field B top at the top of the loop (where a = ir/2). Thus, the 
shape of the loop is given by the equation 



sin a(z) ■ 



B 



top 



B(z) 



(14) 



We can now work out an explicit expression for the field 
line shape z(x) in terms of the function b(z) = B(z)/B(z = 0). 
Let us_ denote the magnetic field strength at the base by 
Z?o = B(z = 0) and the angle between the loop and the verti- 
cal direction at the base by ao(L) = a(z = 0;L). We then have 



b lo JL) = 



Btop(L) 
Bo 



= sinao(L). 



(15) 



Then, the shape of the ascending leg of a loop of length L 
is given by 



x(z) = -b t 



top 



dz' 



b 2 (z')-b 2 op 



L 



top 



dz' 



b 2 (z')-b 



top 

(16) 



The height of the loop, z top , is determined implicitly by the 
condition x(zt op ) = 0: 



top 



dz! 



b 2 (z')- 



-b 2 



(17) 



Ana logo us expressions have bee n obtained by Parker 
(119751) and lBrowning & Priest! (119841) . 

Here are a few analytical examples of the use of this rela- 
tionship. 

1) Example I: 

where Q = z/zo\ zo represents the characteristic magnetic scale 
height of the corona. Then we get b top (L) = 2zo(4zq + L 2 )~ 1 / 2 , 

andz top (L) = (z2 + L2/4) 1 / 2 -z . 

2) Example II: exponential atmosphere, 

b(z) = e~ zlH 

Performing the integration ( TT7| >, we obtain 
—H arctan [e 2(Zlop_z) ' ff - l] 1 ' 2 , and correspondingly, 



(19) 
x(z) = 



L(Zt op ) = 2H arccos [Z?(Z top )] = 2H 



— - ao(Zt p) 



(20) 



Notice that for tall loops with Z top ^> H and b top -C 1, the 
dependence L(Z top ) saturates: L(Z top ) — > L max = ttH. This ex- 
ample illustrates an important point: if the external pressure 
drops off sufficiently steeply, then there is a maximum pro- 
jected length, L max , that a loop in equilibrium can have. This 
implies that if one tries to insert a slender loop with a foot- 
point separation L > L max , then such a loop will not be able to 
attain an equilibrium and will instead grow in height without 
bound, i.e., will tend to open up. 

This fact points to an important feedback: for a given exter- 
nal pressure profile, large loops extend to larger heights, but 
this has an effect of increasing the contribution of these loops 
to the pressure at these large heights (see below), and hence 
may make the pressure profile less steep. 

3.2. Magnetic Energy Density B 2 (z)/8it as a Self-Consistent 
Mean Field 

Our next step is to determine B as a function of z. How- 
ever, the best way to do this is first to express B directly 
in terms of the length of the smallest loop L that reaches 
the given height z, i.e., to find B(L). This can be done di- 
rectly in terms of the orientation-integrated distribution func- 
tion F(L) = J F(L, 6)d0, since, at any given height z, contri- 
butions to the magnetic pressure B 2 /%n come only from those 
loops that extend to this height or higher. Since the depen- 
dence Z top (L) is presumed to be monotonic, then B at a given 
height z will be proportional, roughly speaking, to the integral 
over all loops with lengths L > L(z), where L(z) is the func- 
tion inverse to Z top (L). Thus, naively, we anticipate a result 
that looks something like this: 



B(L)~ A* / F(L)dL 



(21) 



In Appendix A we perform a rigorous analysis and derive 
an exact (within our model) result: 



7rA* - 
db = F(L)dL 



7T A^ 

b(L)=—— I F{L')dL' 



(22) 
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that is, 



Bto P (i) = 7rA* F(L')dL 



(23) 



The condition b(z = 0) = 1, in conjunction with Z top (L = 0) = 
0, gives us the normalization condition for the function F(L): 



F(L')dL' = 



Bo 

7T A* 



(24) 



If all the loops were perpendicular to the disk surface, the nor- 
malization coefficient would be 1/2 (since each loop has two 
footpoints). The extra factor (2/ir) in the above expression 
reflects the fact that small low-lying loops are not perpendic- 
ular to the disk surface, and hence occupy larger horizontal 
projected area on this surface. 

For a given loop distribution function F(L), one can thus 
compute, in principle, the function b(L) that we shall need in 
the next subsection. 

3.3. Self-Consistent Loop Height Z top (L) 

Once L(b top ) is thus determined, we can substitute it into 
equation (T% and thus reduce the whole problem to the fol- 
lowing integral equation for the function U (b) = dz/db: 



Lib) 
2b 



U(b')db' 
^Jb' 2 -b 2 



(25) 



Once this equation is solved, we can integrate U(b) to 
find z(b), thus completing the solution. 

Mathematically, equation d25l l is a linear Volterra integral 
equation of the first kind for the function U (b) in terms of a 
known function L(b). It can be solved exactly. In particular, 
by a simple transformation of variables: t = l-b' 2 , s= l-b 2 , 



G(s) = L(b)/2b,andV(t)-: 
into the Abel equation: 



- U(b') /2b', it can be transformed 



V(t)dt 



■G(s). 



whose solution is 



1 d 

V(s)=- — 

7T dS 




(26) 



(27) 



where we have used G(s = 0) = L(z = 0)/2 = 0. 

Actually, the most useful form of the solution is the first 
equality in equation t27i . By substituting the definitions 
of V(s), G(s) and U(b) into this equation, multiplying by (-2b) 
and integrating with respect to b, we get the following elegant 
final expression for the function z(b): 



z(b) = 



1 f L(b')db' 



T! J \Jb l2 ~b 2 
b 



(28) 



Using equation d23l , we can rewrite this in terms of the func- 
tions F(L) and b(L) as 



Ztop(i) — — — 



A* f L'F(L')dL' 



Bo J ^b 2 (L')-b 2 (L) 



(29) 



3.4. Loop Energy 

One of the main goals of this section is to address the en- 
ergetics of the magnetized corona. Relevant issues include 
the energy distribution of flares as well as the torque on the 
disk due to the coronal magnetic fields. In order to be able to 
address this, we must first determine the energy, £(A), asso- 
ciated with a loop of type A. 

The loop energy is given by the work done by the footpoints 
against magnetic forces as the loop's footpoint separation (i.e., 
the projected length of the loop) is increased from zero to its 
present value L. Since we regard P(z) as isotropic, the energy 
depends only on the length but not on the orientation of the 
loop: 



E(A) = £(L) = / f fv (L')dL' 



(30) 



Here, ff p (L) is the magnetic force on each of the two 
footpoints; it is proportional to the horizontal magnetic 
field Z?hor(z = 0) at the disk surface: 



ML) = ■ 



B-B\ 



C-Dhor 



4tt 



.-=() 



flhorfe = 0) = ^- B hor (z = 0;L) , (31) 

4-7T 



where A^ is the magnetic flux carried by the loop and ah or (z = 
0) = A\I/ /B z (z = 0) is the area of the loop's footspot on the disk 
surface. 2 Thus, we have 



A* 

£(£)=— / B hm (z = 0;L')dL' 



(32) 



Another quantity of interest is the magnetic energy E m!lgn 
contained within the loop: 



« n B 2 (l) 

£magn= / a (I) -jj— dl , 



(33) 



where the integral is taken along the loop from one footspot to 
the other. Using flux conservation, A^>(1) = a(l)B(l) = const, 
this energy can be written simply as 



A* 



B(l)dl. 



(34) 



This expression actually has a very simple physical meaning. 
The integral J B(l)dl is the circulation of the magnetic field 
along the loop; according to the Ampere's law, this is just the 
total surface current flowing around the loop in the perpendic- 
ular direction. Thus, the above expression for ZJmagn is just a 
manifestation of the well-known result that the magnetic en- 
ergy of a current circuit is proportional to the product of the 
magnetic flux enclosed by the circuit and its total current. 

We would like to remark that £(L) can be viewed as a 
magnetic enthalpy, H magn . It includes both the magnetic en- 
ergy E m . dgn stored within the loop and the work W done by the 
loop on the surrounding gas with a fixed (but not necessarily 
uniform) pressure profile: 



(35) 



This is analogous to calculating the amount of heat Q re- 
quired to inflate a hot-air balloon at constant atmospheric 
pressure Pq. Indeed, when the balloon air is heated, the en- 
ergy is expended both to increase the internal energy U of 

2 Note that here, instead of fixing A* and B z (z = 0), we fix and the 
total magnetic field B(z = 0), which includes the horizontal component. 
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the hot air inside the balloon and to perform work against at- 
mospheric pressure (neglecting heat losses from the balloon 
through its skin). Thus, the amount of heat that needs to be 
supplied is equal to the change in balloon's enthalpy H: 



dQ = dH = dU + P dV = 



1 



d(P V)+PodV = -^-P dV , 
7-1 

(36) 

where 7 is the adiabatic index of air. In our case of a magnetic 
loop confined by external pressure, the adiabatic index is 7 = 
2. Therefore, we expect that the total energy that needs to be 
supplied by the footpoint motions to inflate the loop is twice 
the internal magnetic field energy within the loop: 

£ — ^magn W — 2£' ma g n . (37) 

In Appendix B we demonstrate that this is indeed so. 

Note that, in addition to the work W done against the 
external gas pressure as the loop expands and increases its 
cross-section, there is also the work done against the mag- 
netic tension force as the length of the loop is increased and 
the work done on the loop by the external pressure-gradient 
force. However, as long as the loop expands quasi-statically, 
always maintaining its equilibrium shap e, the last two forces 
precisely balance each other (see § I3.lt . and so their corre- 
sponding works cancel out. 

Using the formalism developed in the previous subsections, 
we can now easily calculate the energy associated with a given 
loop. Substituting equation ( fT3l into equation (l32l and using 
expression d23l for B (op (L'), we get 

L L oo 

A* f . A* 2 
£(L) =^J B top (L')dL' = — 

V 

(38) 



7' 



L OO 

J J F(L")dL"dL'. 



The total magnetic energy in the corona is then 

oo oc L oo 

If- „ A* 2 



Eint = — 



F(L)£(L)dL = 



dL / dll \ dL" F(L)F(L") 



o 



A* 2 



o 



J dL J dL' L' F(L)F(L'). 



(39) 



It is instructive to consider a case in which F(L) has a power 
law tail, F ~ L~ a , truncated at some large L max 3> H. Then, as 
can be seen from equation d38l . the energy of the largest loop, 
£(L max ), is almost independent of L max for a > 2, but starts 
to grow as L 2 ^™ for a < 2. In a similar manner, from equa- 
tion ( |39l it follows that the total coronal energy becomes dom- 
inated by the large-L tail if F(L) drops off sufficiently slowly, 
i.e., a < 3/2. In this case the total energy scales with L max as 
E to t oc L^ 2 " and may become much larger than the potential 
magnetic field (whose characteristic magnetic scale height is 
of order H). Physically, we expect L max to be at most about 
the local disk radius, r, so that E tot is bounded. 

Finally, let us consider the angular momentum transfer by 
the coronal magnetic field. The torque due to a single loop is 
given by 

A* AyAr AVB 

#hor — ; — = 5 b top (L) L sin 20, 

L 07T 

(40) 

and hence the total torque per unit disk area is 

G , = _A*Bo j [jLdOFmfibtopmLsinie. (41) 



AG(L,0): 



4tt 



_ Let is again consider the truncated power-law example, 
F ~ L~ a( - e K For a fixed degree of anisotropy, e.g., a fixed 
the characteristic angular scale | sin 8 m \ n | at which the function 
a(0) has a minimum a m i n , the torque becomes dominated by 
large loops if a m \ n < 3/2, similar to E tot [where we used equa- 
tion ( 1231 . It then scales as L^~ 2 " and thus may become by a 
factor (L max /H) 3 ~ 2a ^> 1 larger than the usual torque exerted 
directly by MRJ turbulence within the disk, assuming that the 
magnetic field at the base of the corona is comparable to that 
in the disk. In reality, however, a decrease in a m \ n may come 
hand-in-hand with an increase in the degr ee of anisotropy of 
the distribution function F(L, 9) (see § 15. 3t . manifested as a si- 
multaneous decrease in | sin6* m j n |. If this is the case, the torque 
amplification will be not as strong. 

4. THE LOOP KINETIC EQUATION 

In this section we discuss how to calculate the loop distri- 
bution function. In particular, we construct the loop kinetic 
equation that governs the evolution of this function. 

4. 1 . Physical Assumptions of our Model 

In order to build a quantitative model of the magnetic field 
in the corona, we need to make some specific assumptions 
about the most important physical processes that govern the 
life of individual coronal loops, including their interactions 
with each other. These assumptions are the main building 
blocks of our model; we shall discuss them in this section. 

1) The Alfven velocity in the corona is much faster than 
both the disk's rotational velocity at the given radius and the 
thermal velocity of the coronal gas; therefore, the corona is 
considered to be in a slowly-evolving force-free magneto- 
static equilibrium at all times and almost everywhere (except 
for rapid rearrangements due to reconnection events, see be- 
low). 

2) The disk is geometrically thin, with the gas scale-height 
much smaller than the distance from the central object, H <C 
r. This gives us an important small parameter that can be 
used in the analysis. For example, this assumption gives us 
an "inertial range" of spatial scales much smaller than r but 
much larger than H. This enables us to perform an analysis 
that is local in r. In particular, this means that we can neglect 
geometrical effects resulting from cylindrical geometry when 
considering the flux-loop expansion process. Note, however, 
that the validity of the thin-disk assumption is questionable 
close to a black hole accreting near its Eddington limit. 

3) The disk is differentially-rotating (e.g., keplerian). As 
a result, coronal loops with radially-separated footpoints are 
subject to continuous stretching in the toroidal direction. This 
generates toroidal magnetic field whose pressure inflates the 
loops and ultimately leads to the creation of a vertically- 
extended corona (see below). 

4) At any given time, the shape and the overall height of 
each loop are determined by the magnetostatic equilibrium of 
the loop as if it were confined by a stratified atmosphere with 
certain external isotropic pressure P(z). In turn, this pressure 
represents the effective magnetic pressure of all the neighbor- 
ing loops, and we shall denote it as P = B 2 {z)/%tt. This equi- 
librium shape is maintained at all times, since it adjusts on the 
Alfven time scale, which is assumed to be much faster than 
the disk rotation. 

5) The disk is turbulent due to the usual internal MRJ (as 
opposed to coronal that may act simultaneously). The char- 
acteristic spatial scale of the turbulence is H, and the char- 
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acteristic time scale is f^. The important effects of the disk 
turbulence on the corona are: 

5a) Flux Emergence plays a very important role in the so- 
lar corona and, by analo gy, is also b e lieved to be important 
in the case of the disk Gale ev et al.l (|T979). We generally 
expect the emerging magnetic loops to be relatively small 
in size (of order the disk thickness H) and to have typical 

1/2 

magnetic fields of order a ss ' B ea , where the dimensionless pa- 
rameter a ss ~ 0.01 -0.1 is the Shakura & Svun vaevl (1 19731) 
viscosity coefficient and B eq is the field strength that corre- 
sponds to equipartition with the gas pressure inside the disk. 
In addition, numerical studies of MRI turbulence show that 
the toroidal field in the disk tends to be larger than the radial 
field by a factor of 5-10. Thus, flux emergence is expected 
to be anisotropic, with newly emergent loops elongated in the 
toroidal direction by a factor of a few. 

5b) Again, similar to the Sun, the disk turbulence leads to 
a two-dimensional random walk of the coronal loops' foot- 
points on the disk surface. We expect this random walk to 
be characterized by spatial and temporal scales of the order 
of H and respectively. However, similar to the process of 
flux emergence discussed above, the random walk, in general, 
may be anisotropic, with characteristic steps in the azimuthal 
direction being somewhat larger than in the radial direction. 

6) Reconnection: In our model, two loops may reconnect 
with each other, forming two new loops (see Fig . O ■ Thus, re- 
connection represents a binary interaction between individual 
magnetic structures, analogous to binary collisions between 
particles in a gas. 

We shall assume that, once triggered, a reconnection event 
(a flare) happens very quickly, essentially instantaneously on 
the orbital time scale. This assumption can be justified by not- 
ing that the corona is assumed to be a very low-density, and 
hence collisionless, environment. Therefore, reconnection 
there proceeds in the Petschek-like fast collisionless regime, 
enabled by anomalous resistivity or by the two-fluid (Hall- 
MHD) effects. The characteristic reconnection time-scale is 
then only by a factor of 10-100 slower than the Alfven cross- 
ing time ta- Thus, since we assume that Va ^S> Vk, it is rea- 
sonable to expect that the typical duration of coronal recon- 
nection events may still be fairly short compared with the or- 
bital time-scale VT l . Then, to the extent that f2 _1 is the main 
dynamical time-scale in our problem, characterizing differ- 
ential rotation, flux emergence, and turbulent random walk, 
we can, for the purposes of our study, regard reconnection 
between loops, once triggered, as being essentially instanta- 
neous. Thus, we arrive at a picture in which magnetic loops 
evolve slowly {i.e., on the orbital time scale), but from time to 
time they suddenly and instantaneously reconnect. This pic- 
ture is similar to the observed behavior of solar coronal loops, 
where the characteristic reconnection (or flare) time is typi- 
cally much shorter than the typical loop lifetimes. Thus, from 
the standpoint of viewing the corona as an ensemble of many 
loops, reconnection events can be regarded as relatively infre- 
quent binary collisions between loops, analogous to the binary 
collisions between particles in Boltzmann's gas. An impor- 
tant corollary from this is that the footpoints of the loops do 
not have time to move significantly during the reconnection 
event. As we shall see in § 14.21 this will give us the rules that 
determine the footpoint separations of newly-formed loops. 
We shall also assume that these newly-formed loops quickly 
assume their equilibrium shapes (see above). 



reconnection point 




FIG. 5. — Reconnection between two loops as a binary process. 



4.2. The Loop Kinetic Equation 

Based on the above assumptions, we shall now work out an 
evolution equation for the loop distribution function F . We 
shall call this the Loop Kinetic Equation (LKE). This equa- 
tion should have the following terms, reflecting the relevant 
physical processes: 

1) Flux emergence/submergence acts as a source/sink of 
new coronal loops. It can be modeled by a source term, 
S{A), that describes the rate at which the loops emerge into 
the corona, their characteristic sizes, magnetic field strengths, 
etc., (or, in a more elaborate model, by specifying the distri- 
butions of these quantities). Specifically, one can add loops 
at some characteristic "injection scale", somewhat larger than 
the disk thickness H, and remove very small loops, say, of 
size H or less, as it was done in the model by Hugh eset alj 
(120031) . In addition, in the ADC case, we expect flux emer- 
gence to be anisotropic, with the emerging loops being by 
a factor of a few longer in the toroidal direction than in the 
radial direction, a s indicated by numerical simulations {e.g. 
iHirose et al.ll2006h . 

An alternative way to take flux emergence into account is 
via the boundary conditions for F at small scales (of order H). 
This is the view adopted in our model. This choice is justified 
by arguing that the population of smallest loops comprising 
the "magnetic carpet" is predominantly determined by a de- 
tailed balance that is quickly established with the magnetic 
fields in the disk itself. This process turns over (operates) 
very quickly and hence the distribution of the very small loops 
is basically independent of what happens in the larger-scale 
corona. 

2) Random footpoint motions due to the disk turbulence. 
Since we expect the characteristic steps of this random walk 
to be relatively small (of order H, see § 14.11 ) compared with 
the sizes of most loops under consideration, we can employ a 
Fokker-Planck-like approach to this process. This results in a 
dif fusion operator with the diffu sion coefficient of the order of 
thelShakura & Svun vaevl (I 19731) a- viscosity coefficient: D ~ 
a ss c s H ~ otssQH 2 . In general, however, this diffusion may 
be anisotropic, with D being a tensor {e.g., a diagonal tensor 
with Dyy > D rr ). We expect the effect of the random walk to 
be relatively unimportant for large loops, L> H. 

3) Keplerian differential rotation leads to a secular evolu- 
tion of Ay, 



Here, SI = il^(r) can be regarded as constant because we con- 
sider spatial scales that are small compared with the disk ra- 
dius, Ar <C r. In the Loop Kinetic Equation this process is 
described by an advection term, (3/2) £1 Ar{dF / dAy). 

4) Coronal MRI backreaction term ^backreaction : m a 
geostrophic approximation, this is obtained by balancing the 
magnetic force on the footpoints (per unit area) with the Cori- 
olis force (also per unit area) due to the rotation of the loop it 
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induces: 



2 gr . Bhor = 2S[0x £v backreaction ] , (43) 

47T 

in which Vbacfaeaction represents the departure from the keple- 
rian rotation velocity. 

5) Interaction of two loops by reconnecting with each other, 
yielding two new loops. This process is described by a binary 
collision term F lec . This is the most non-trivial term and we 
devote the entire next subsection (§ |4.3t to a detailed discus- 
sion of it. 

At the end, we arrive at the following general form of the 
kinetic equation for the loop distribution function: 

dF ( <9 2 <9 2 

-(Ar, Ay, t) = S(Ar, Ay) + D rr 1 r lr ^+D yy I F 



dt 



dAr 



dAy 1 



3 dF 

+ - f2 Ar — — + Fbackreaction + F rec (44) 

2 oAy 

The simplest meaningful case of this equation is when one 
neglects the source, diffusion, and feedback terms and looks 
for a steady state that is produced by the balance between 
Keplerian shear and reconnection: 



2 dAy 



(45) 



will thus be analogous to, and can 
extension of, the previous work by 



Our approach 
be regarded as an 

iTout & Prin gle ( 1996), whose main goal was to study the for- 
mation of large magnetic structures via the reconnective "in- 
verse cascade" in the corona. Following them, we also repre- 
sent the coronal magnetic field by an ensemble of flux loops 
described by a distribution function. However, our model is 
more general and uses more realistic p hysics. We take int o 
account a number of effects ignored in lTout & Pringld dl996l) . 
such as inflation of the loops as they are stretched by the Ke- 
plerian differential rotation. Also, in their model reconnec- 
tion was taking place only at the disk surface and thus one of 
two newly reconnected loops was vanishingly small and was 
assumed to just disappear; as a result, the reconnection pro- 
cess did not conserve the number of loops. We, on the other 
hand assume that reconnection occurs higher in the corona, 
and hence two new loops form and the loop numb er is con- 
served (similar to the model by Hughe set alj (120031) ). 

4.3. Reconnection Described as a Collision Integral 

Two loops A and B may interact by reconnecting with each 
other and forming t wo new loops C and V as a result (see 
Fig- IS). Following ITout & Pringld (119961) . we shall describe 
this process by a nonlinear binary-collision integral, similar 
to the collision integral in the Boltzmann kinetic equation. In 
reality, of course, interaction between loops is more compli- 
cated and so such a description is oversimplified. Moreover, 
magnetic loops fill up the entire coronal space, and so they 
resemble more a non-ideal liquid rather than an almost ideal 
rarefied gas with infrequent binary encounters. Nevertheless, 
we believe that this binary-collision representation of recon- 
nection can lead to some valuable physical insight into the 
complicated dynamics of the coronal magnetic field. 

The Boltzmann collision integral can be split into two: the 
source term and the sink term. The sink term F C0 \\-(A) de- 
scribes the rate of reduction in the number of loops of a given 
type A due to reconnection between these loops and all other 
loops. The source term F co \i+(A) describes the rate of increase 



in the number of loops of type A when they are a product of 
reconnection of other-type loops. By "type" we here mean a 
set of loops with the same values of their primary parameters 
(L, 9) or ( Ar, Ay). Thus, each of these terms is a quadratic in- 
tegral operator, with a kernel that depends both on the types of 
the two loops and also on their relative position (see below). 
Thus, we can write the reconnection term schematically as 



I*rcc(.A) — F rsc —(A) 



where 



F KC .-(A)-- 



M)=: 



dB Q AB F( A)F(B), 



(46) 



(47) 



dCdVQ cv ^ A F(C)F(V). (48) 



Here, dB = dArgdAyg, etc., The kernels Qab in the sink 
term and Qcv^A m the source term are related via 



Qab = ^ / llL ' Qab- c ■ 



(49) 



Using this relationship, the two terms can be combined as 



dCdVF(V) Qcv->aF(C)-Q A v^cF(A) 

(50)" 

In order to go from this general expression to a specific 
operational procedure, we need to formulate the rules that 
govern the reconnection process. Indeed, the two new loops 
C and T> formed as products of reconnection between two 
loops A and B cannot be arbitrary and certain selection rules 
must be applied. More specifically, from the assumption that 
reconnection is instantaneous on the orbital time scale, it fol- 
lows that the footpoints of the reconnecting loops do not move 
significantly during the reconnection event. Only the way they 
are connected to each other changes. Therefore, the primary 
parameters (i.e., footpoint separations) of loops C and T> are 
uniquely determined by the footpoint positions of the two in- 
coming loops A and B: 

(51) 
(52) 



Ar c = r c - - r c + = r A - - r B+ , 
Ar D = r D - - r D + = r B - - r A+ , 



and similarly for Aye, Ayo- Here, A+, A_ are the positive and 
negative polarity footpoints of the loop A, etc. Mathemati- 
cally, these rules play a role similar to the momentum and en- 
ergy conservation conditions for particle collisions in kinetic 
theory of gases; they enter as (5-functions in the interaction in- 
tegral in our loop kinetic equation. Thus, one can easily see 
that the parameters of the new loops depend not only on the 
parameters of the old loops, but also on the positions of A 
and B relative to each other (see below). 

The kernel Qab m equation (l47l i is the probability rate {i.e., 
probability per unit time) that two loops of types A and B will 
come together and reconnect. Thus, Qab describes the rate of 
reconnection events (the number of such events per unit time). 
This should not be confused with the concept of "reconnec- 
tion rate", a widely-used term in reconnection research with a 
completely different meaning. 

Now let us discuss on which physical parameters Qab 
should generally depend. Whereas Tout & Pringle (1996) 
just took Qab - const, we want to develop a more realis- 
tic and more sophisticated model, taking into account several 
important factors. First, notice that Qab has dimensions of 
[cm 2 /sec]. Based on dimensional arguments, it should then 
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be proportional to the characteristic rate at which the coro- 
nal magnetic field is reconfigured. In addition, it should also 
reflect the fact that larger loops have larger "interaction cross- 
section" (see below), and thus should be roughly proportional 
to the squares of loop sizes. 

Let us first address the characteristic reconfiguration time- 
scale. The only fundamental dynamical time-scale in the 
corona, i.e., the time-scale on which the corona, seen as an 
ensemble of elementary coronal structures, rearranges itself, 
is the orbital time, fi" 1 (or the inverse of the shear rate, 3/257, 
which is not independent). This means that, if we represent 
the evolution by a sequence of discrete steps, each step rep- 
resenting a noticeable change in the relative position and/or 
orientation of the coronal elements, then the most appropriate 
choice for duration of these steps is of order fi . Therefore, 
in general, Q should scale with ft. Next, if we follow a given 
magnetic element, at each new step there will be a certain 
probability n < 1 that the resulting new magnetic configura- 
tion around this element becomes favorable for reconnection 
of this element with another. We shall treat k as a constant 
number, independent of the loops primary parameters. Thus, 
the overall rate at which the loops are disrupted through re- 
connection with other elements should be proportional to kQ: 
Qab = Ktt<j AB , (53) 
where we introduced the "reconnection cross-section" cfab- 

The cross-section o A b should in some way scale with the 
loop sizes. It involves contributions from all possible relative 
positionings of the two interacting loops for which the two 
loops "effectively intersect". We shall describe this relative 
positioning by two impact parameters, bu and b±, defined as 
the offsets between the centers of the two loops in the direc- 
tion parallel and perpendicular to loop A, respectively. We 
shall assume that once the impact parameters are in a range 
such that the two loops "effectively intersect", the probabil- 
ity that these loops will reconnect is constant, independent 
of their positions or their parameters. Furthermore, for sim- 
plicity we shall assume b^ and b± to be uniformly-distributed 
independent random variables. Thus, uab is just equal to the 
"interaction area" in the (b|| ,b±) space that corresponds to an 
"effective intersection" of the given two loops: 



cab = / da AB = 



db\\ db \ 



(54) 



But what do we mean by "effectively intersecting"? If 
the two loops are approximated by their central lines (one- 
dimensional objects), then the set of values (b|| , b±) for which 
they intersect is also one-dimensional (a line segment), and 
thus has measure zero in the two-dimensional (b^,b±) space. 
In other words, the probability that two randomly-drawn lines 
intersect is zero in the three-dimensional space. Therefore, 
to get a meaningful result, we need to take into account fi- 
nite thicknesses of the loops. In particular, we shall say that 
two loops "effectively intersect" when the closest distance be- 
tween their central lines is less than a certain fraction of the 
combined loop thicknesses at the intersection height. Oper- 
ationally, for a given value of bu, say, we can introduce the 
y-cross-section (JAB,±(b\\) as the spread in the values of b± 
which result in an effective intersection of two given loops A 
and B. We can then write 

<Jab = J <T A B,±(b\\)db\\ . (55) 

where the integral is taken over the range of impact parame- 
ters bu for which an intersection between loops A and B is at 
all possible. 



Correspondingly, the sink-term part of the reconnection in- 
tegral can be written as 

F KC AA) = -nn J J dBdb\\a A B,i_{b\\)F(A)F{B). (56) 

Now let us consider the source term (t48t . Employing the 
arguments given earlier in this subsection, we introduce the 
cross-section for two given loops C and T> to reconnect giving 
a loop A as a result: Qcv^A = nil (Jcd^a- We can then write 



dC dV (Jcv^a F(C) FiV) . (57) 



The V integral in this equation is taken over all the loops 
that can yield a loop of type A as a result of reconnection 
with a loop of a given type C. Note that, whereas one did not 
need to know reconnection product loops to compute the sink 
term ([56}, in order to calculate the source term, this knowl- 
edge is in fact necessary. It is contained in our "reconnection 
rules", such as those given by equations (T5TT> — <T52l>. For def- 
initeness, let us consider the case when the resulting loop A 
starts from the "+" footpoint of loop C and ends at the "-" 
footpoint of loop T>. Then, for given C and A, there is a well- 
defined range of impact parameters b\\ for which one can find 
one (or sometimes two) loop T> = T> A -c(b^) that intersects 
loop C and ends at the end footpoint of loop A. Thus, in 
principle, for a given b\\ one can formulate the rules that re- 
late the primary parameters Aro and Ay a to those of loops C 
and A. We shall denote these relationships by Ar^~ c (/7 r ) and 
Ay^~ c (fo||). In general, there may be one or two such solu- 
tions. 

Recalling now that flux tubes have a finite thickness, we 
have, by analogy with equation i 



<Jcd~,a =2 J db\\ c7 C D.Ab\\) S[V-V A - C (b\\)] . 



(58) 



where 5[V -V A . c {b\\)] = 5[Ar D - Ar^ib^)} x 5[Ay D - 
Ay£," c (b||)]. The factor 2 accounts for the fact that in the pre- 
ceding paragraph we considered only one half of all possible 
configurations, requiring the starting footpoint of loop A to be 
the starting footpoint of loop C. For each such configuration 
there will also be an identical contribution from interchanging 
loops C and D. 

Substituting this cross-section into our expression d57l i 
for F co \i + (A), and using the (5-function to integrate over dT> = 
dAr^dAyo, we get 



F con , + (A) = Kfl J J dCdb {l a C D^c.Ab\{)F{C)F[V A . c {b\\)}. 

(59) 

Combining equations d56l ) and d59l . we can write: 
F rec (A) = tfi / dC dbu F(C) 



F[V A - C {b\\)]acD A ^Ab\\)-F{A)aAC.±.{b\\)m 



Effect of Finite Loop Thickness 

For a given value of bu, one first finds the coordinates 
[r(b\\),y(b\\),z(b\\)] of the point where the central lines of the 
two loops would intersect. After that, one calculates the _L- 
extent of each of the two loops at this point, which hence gives 
one (7 AB \{b\\). 
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When doing this, one should take into account the follow- 
ing important effect. If the loops have no internal twist, as we 
assume here, the longitudinal magnetic field is approximately 
constant across a loop and is roughly equal to the characteris- 
tic magnetic field B(z) at a given height z. Then, by flux con- 
servation, the cross-section of the loop (normal to its central 
line) varies along its length, and, in particular, may increase 
greatly at large heights if B(z) drops off rapidly. More pre- 
cisely, approximating a loop's cross-section at a given height z 
as a circle of some radius d(z), we can estimate this radius 
as d(z) ~ [A*/7rB(z)]'/ 2 = d Q [B(z)/B Q Y l l 2 . As long as the 
loop is slender, d ^L, this should be a good approximation. 
Thus, if the intersection point lies high above the disk, so that 
B(z) <S Bo, the characteristic loop thickness is much larger 
than near the disk surface. Consequently, the reconnection 
cross-section is increased, which has important implications 
for the "inverse cascade" of magnetic loops in the corona. It 
is also interesting to note that this process is controlled by the 
self-consistent field B{z). 

[Note that, in the case of twisted loops this effect is not 
as profound and a more accurate approximation is probably 
given by d(z) — do = const (see § I6.21 >.1 

Thus, typically we expect <tab,±. to scale as 

(61) 

For example, if B is the smaller of the two loops, we expect 
the typical interaction height to be of the order of this loop's 
height Zb- Correspondingly, we expect ctab,± ~ dob~ l l 2 {ZB), 
where b(z) = B(z)/Bq. This estimate for <TAB,±(b\\) will be 
roughly valid for almost the entire allowed range of bu, which 
is of order 4 Ar^. 

To sum up, larger loops have larger cross-sections for re- 
connection, for two complementary reasons. First, the cross- 
section is enhanced because larger loops have larger range 
of impact parameters (in the radial direction, say) for which 
intersection of their central lines is possible. Second, larger 
loops extend to, and may interact at, larger heights, where the 
mean magnetic field is weaker and hence where the loops be- 
come fatter. They thus have a greater chance of overlapping 
with each other. As a result, the reconnection cross-section 
scales with the size Lb of the smaller of the two loops as 

a A B~LBdob-' /2 (Z B ), (62) 

[for Lb ^> d(Zg)]. We thus see that the function B(z) affects 
the evolution^>f the loop distribution function. Since, accord- 
ing to § 13.21 B(z) is itself determined by the distribution func- 
tion, this means that determining these two functions together 
in a self-consistent way requires an iterative procedure. 

What is important here, is that larger loops have a tendency 
to reconnect with each other quickly, probably leading to a 
rapid "inverse cascade" to even larger loops. In addition, the 
fact that reconnection with large loops cannot be neglected 
suggests that a Fokker-Planck-like approximation to the re- 
connection term will not work. This is because reconnection 
events that lead to large changes in loop parameters are im- 
portant and so our collision integral cannot be described by a 
differential diffusion-like operator. 

5. NUMERICAL SOLUTION 
5.1. Numerical Setup 



We solve the LKE numerically and obtain a steady state 
solution. For simplicity we leave out the turbulent diffusion 
term, the coronal backreaction term, and the source term (flux 
emergence is then treated via the boundary condition at small 
scale, see § 15.2b . Thus, we include only the two processes 
which we believe are the most important: Keplerian shear and 
reconnection between loops. Correspondingly, we aim here at 
investigating the effect of the relative importance of these two 
processes on the steady-state loop distribution function. 

In the numerical implementation we work in the (L, 8) pa- 
rameter space where L = (Ar 2 + Ay 2 ) 1 ' 2 is the distance be- 
tween the footpoints and 6 is the angle the vector (Ar, Ay) 
makes with the toroidal direction, measured clockwise: 
tan# = Ar/Ay. We use a grid that is uniform in 9 (between 
and 2tt) but logarithmic in L. The L grid spans from L^ = 1 
to some L max ^> 1 (usually we take L max to be 10 or 20) in 
length units such that H ~ 1 . 

The advection term resulting from Keplerian shear is very 
easy to implement, we just use one-sided derivatives. 

The reconnection collision integral is obviously more com- 
plicated and we devote the rest of this subsection to our nu- 
merical implementation of it. At each time step we go over 
all possible pairs of loops A and B. Furthermore, for each 
given pair we go over all possible reconnecting configura- 
tions, distinguished by the impact parameter bn, defined as 
the displacements between the centers of the two loops in the 
direction along loop A (see § I4.31 l. Thus, at each timestep, 
we are performing a five-dimensional (5D) integration, which 
makes increasing resolution extremely numerically costly. 
We found, however, that numerical convergence is very good 
and relatively modest resolution suffices. To study numerical 
convergence, we performed calculations with Nl = 20 and 40 
points in logL, Ng = 20, 40, 60, and 80 points in 6, and 30, 
50, or 100 points in b\\. We found the resulting F{L,8) to be 
essentially unchanged, although small values of n required a 
higher 0-resolution for convergence (see below). 

For simplicity, in our treatment of reconnection, we assume 
the loops to be semi-circular in s hape, instead of using equi- 
librium shapes discussed in § 13.11 This simplification has two 
benefits. First, since the shapes are described by relatively 
simple analytical expressions, we can derive explicit analyti- 
cal relationships expressing the parameters of the two product 
loops in terms of the parameters of the interacting loops and 
the impact parameter (we call these relationships the "recon- 
nection rules"). Having such expressions in an explicit form 
greatly simplifies the numerical procedures. 

The second advantage of the semicircular approximation 
stems from the observation that, because all the loops have the 
same shape, the reconnection rules essentially depend only 
on the ratio of loop sizes, whereas their dependence on the 
absolute loop size is a trivial rescaling. Likewise, the recon- 
nection process is, by itself, isotropic, i.e., does not depend 
on the absolute orientation of the loops, only on the angle 
between them. This enables us to to reduce the analysis of 
reconnection between two given loops A and B to consider- 
ing a template that corresponds to the given angle between 
the loops and the given length ratio Lb /La- That is, we can 
analyze reconnection between any two loops in a rotated and 
rescaled system of coordinates, (x' ,y'), in which one of the 
loops (loop A for definiteness) is in the positive x'-direction, 
i.e., has 9' A = tt/2 and has unit length, L' A = 1. In practice, 
we first (even before we start the evolution of LKE), create 
a lookup table describing the reconnection rules for interac- 
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tion of this loop A' with all other loops £>'. The lookup table 
is 3D, two of the coordinates being 0' B and L' B , and the third 
coordinate being the rescaled impact parameter b'« ; the lat- 
ter lies within the range from -(1 + Ax' B )/2 to (1 + Ax' B )/2, 
where Ax B = L' B sm9' B . For each given 9' B , L B and b',,, it is 
just a matter of simple algebra to figure out the perpendicular 
displacement b' ± corresponding to the intersection between 
the semi-circular loops. Simultaneously one finds the 3D po- 
sition of the intersection point between loops A' and B 1 for 
given bj| (including the height of the reconnection point). One 
can then readily deduce the parameters of the two product 
loops C and V . It is easy to see that the problem reduces to 
quadratic equations and hence for each 9' B , L B , and b'n there 
may be zero, one, or two solutions. 

During the time evolution we use this table as follows. For 
each given pair of loops A and B we find their corresponding 
template pair by rotating them by 9a~tt/2 and rescaling the 
loop sizes by La- We also rescale the impact parameter by La, 
i.e., bj| = b\\ I La- We then use the template table to find the two 
template product loops and we transform them back by multi- 
plying by La and rotating by w/2- 9a to find the actual product 
loops C and T>. Having a grid uniform in 9 and logL makes 
this procedure especially convenient and straightforward. 

Once the product loops are found, one also needs to figure 
out the reconnection cross-sections corresponding to given A, 
B, and b'^ . In accordance with the above discussion, daAB(b^) 
is proportional to db\\ = La (I + Ax' B )/Nb) and to the loop 
combined deprojected thickness db' ± at the intersection point. 
Apart from simple geometrical projection factors, db' ± is pro- 
portional to d(z) ~ dolBo/Biz)] 1 / 2 , where z is the height of the 
reconnection point above the disk. This factor enhances the 
reconnection probability for large loops intersecting at large 
heights. Since B(z)/Bo itself depends on F(L,9), this proce- 
dure requires iteration. Once the reconnection cross-section 
is found, one can proceed to evolve the number of 
incoming (A and B) and product (C and T>) loops. Namely, 
F{A)dA and F(B)dB are reduced by 

F(A)F(B)dAdBL A db\ l KilaMb'^At , (63) 

and F(C)dC and F(D)dT> are increased by the same amount 
at each timestep Af . 

Finally, we varied the initial conditions and found that our 
resulting steady state solutions are insensitive to them. 

5.2. Boundary Conditions 

For various physical reasons, we expect our model to break 
down at both small scales and large scales. We thus need to 
discuss how to prescribe the boundary conditions for the LKE 
at both of these scales. 

At small scales, L ~ H, the model is expected to become 
invalid because the magnetic field is no longer force-free, as 
plasma pressure starts to become dynamically important at 
small heights. In addition, the characteristic thickness of mag- 
netic structures near the disk surface is expected to be of or- 
der H and hence small loops with L ~ H are not going to 
remain slender as our model assumes. Thus we need to un- 
derstand what is a plausible way to describe these small loops. 
At small scales the dynamics of magnetic loops is strongly af- 
fected by rapid flux-exchange {e.g. flux-emergence) processes 
with the turbulent disk. One can therefore argue that the over- 
all distribution of the smallest loops is largely determined by 
the detailed balance equilibrium that is rapidly established be- 
tween smallest loops and the turbulent disk. Characterizing 



this intense interaction with the disk turbulence and calculat- 
ing the small-scale distribution function is completely beyond 
our present model. It would require thorough understand- 
ing of 3D MHD turbulence in stratified disks with simultane- 
ous a ctions of the MRI and Parker instability dTout & P ringle 
1992); most likely, this task will be accomplished using so- 
phisticated numerical simulations. However, what matters for 
our model here is that we expect that the number of small 
loops be essentially insensitive to what happens to larger coro- 
nal loops. This means that we can mimic the disk-corona in- 
teraction by setting-up a Dirichlet-type boundary condition at 
small scales, i.e., by prescribing the distribution function at 
some small-cutoff scale L mm ~ H: 

F(L = L mm ,9) = F l (9). (64) 

(In our model we set L mm = 1.) Notice that F\(9) is in gen- 
eral expected not to be isotropic because emerging loops 
will be preferentially azimuthally elongated (| sin#| <C 1). In 
the present simulations, however, we take it to be isotropic, 
Fi(9) = const. 

Our model also breaks down at large scales. In particular, 
when a loop's size becomes comparable to the disk radius r, 
the local cartesian geometry adopted here is no longer appli- 
cable, and toroidal effects become important. This leads to a 
much faster expansion of stretched loops, resulting in a com- 
plete opening of the field, as discussed in § 12.21 In the future, 
we plan to implement a physically realistic way of treating the 
opening of large loops by incorporating open field lines into 
our model (see § |6.31 >. 

In the present model, however, we just introduce some large 
cut-off scale, L max , and we set boundary conditions at this 
scale. We experimented with two types of boundary condi- 
tions: in one case, we set F(L max ) = 0, i.e., we simply remove 
all loops that reach the cut-off scale. In the other case, we just 
limit the growth of loops beyond L max ; i.e., as soon as a loop's 
length exceeds L max , we reset it back to L raax . This leads to a 
gradual pile-up of loops near L max , but does not affect the loop 
distribution on smaller scales; in particular, it does not change 
the power law tail for sizes just slightly smaller than L max . 

5.3. Results 

Distribution Function 

We performed a series of calculations with several differ- 
ent values of k (k = oo, 1, 0.3, 0.1, 0.03, 0.02, 0.01, 0.005, 
and 0.003). Some of the resulting steady state loop distri- 
bution functions for our fiducial resolution Nl = 40, N$ = 80 
and for the large-scale boundary condition F(L m!lx ,9) = 0, 
Lmax = 10, are presented in Figures [6H5] In particular, Fig- 
ures [6] and [7] show (in log-log coordinates) F as a function 
of L for purely radial loops (9 = n/2) and purely toroidal loops 
(9 = 0), respectively. In the case k = oo (keplerian shear turned 
off), the distribution function is isotropic. This is of course 
expected, since loop-loop reconnection — the only process 
determining the distribution function in this case — by itself 
is independent of the absolute orientation of the reconnect- 
ing loops. As the frequency of reconnection events relative to 
shear, quantified by k, is decreased, a given loop (especially if 
it has a large |Ar|) experiences, on average, larger stretching 
in the toroidal direction by the shear before it undergoes re- 
connection with another loop. As a result, the loops become 
predominantly azimuthal and the loop distribution function 
becomes more and more anisotropic: F(9 = w/2,L) steepens, 
whereas F(9 = Q,L) becomes shallower with decreasing k. 
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Because our problem lacks a preferred length-scale be- 
tween L m ; n and L max , we find that, along each ray 9 = const, 
F{L) is well described by a power law with the orientation- 
angle-dependent exponent: 

F(L,6) ~ L- a ^ 0) . (65) 

Figure [8] presents the function a K (9) for k = 1, 0.3, 0.1, 
and 0.03. 

Finally, in order to estimate the total magnetic energy in the 
corona and the total magnetic torque (see § 13 .41 . one needs to 
know the (^-integrated distribution function, 

2tt 



F(L)= / F(L,9)dd. 



(66) 



In Figure we plot F{L) for n = 1.0, 0.3, 0.1, 0.03, 0.01, 
and 0.005 in log-log coordinates. In general, of course, one 
cannot expect a (^-integral of exponents L~ a ^ to be itself a 
power law of L. However, we find that the integral is strongly 
dominated by the range of values 9_ corresponding to the ridge 
in F(L, 9), and so a power law F(L) ~ L _Q ( K ) is actually a 
reasonably good approximation, especially for relatively large 
values of k. Because the power-law fit is not ideal, there is 
some degree of uncertainty in determining the value of a; we 
estimate the characteristic error to be < 0.1. 



The dependence of the power-law exponent a on the re- 
connection parameter k is plotted in Figure [10] We find that 
when reconnection is strong compared with keplerian shear 
(k > 0.2), the power-law exponent stays close to 3.1, inde- 
pendent of k. This is because, as long as the keplerian shear 
can be neglected, reconnection is the only term on the right- 
hand side of the LKE. In this case, k cannot affect the steady 
state solution, it can only regulate how fast this solution is 
achieved. As k is decreased, however, the keplerian shear 
term becomes important and a(n) starts to decrease. We find 
that its overall behavior can be approximated by o.(k) ~ 3. 1 = 
const for k > 0.2 and <5(k) ~ 3.55 + (1/3) In n for k < 0.2. 
Correspondingly, we expect a(n) to crosses the critical val- 
ues 2 and 3/2 (see §[331 at k 2 ~ 0.01 and k 3/2 ~ 0.002, re- 
spectively. As a reminder, a < 3/2 means that the total mag- 
netic energy of the corona is dominated by the largest loops. 
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FIG. 9. — Orientation-averaged loop distribution function F(L) for k = 1.0, 
0.3, 0.1, 0.03, 0.01, and 0.005. 



Finally, although most of our runs were done with L max = 
10, we also conducted some runs with L max = 20. We found 
the distribution function to be essentially the same; the power- 
law just extended further in the L max = 20 case but the slope 
and the normalization were unchanged. 

Energetics and Torque 

Using the computed loop distribution functions, we calcu- 
late the magnetic energy density as a function of height above 
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FIG. 10. — Power-law exponent a for the orientation-averaged loop distri- 
bution function F(L) ~ L~ a , as a function of re. The dotted horizontal lines at 
o = 2 and a = 3/2 correspond, respectively, to the critical values at which the 
magnetic energy of the largest loop, E(L milx ), and the total magnetic energy 
in the coro na, £ t ot, start to become dominated by the large-scale cut-off L m;lx 
(see §|H. 



the disk, B 2 (z)/8ir, and the energy of a loop as a function of 
its length, £(L), for several different values of k (see Figs.fTTI 
and[T2l. We see that these functions remain essentially inde- 
pendent of k as long as it is large enough, k > 0.1), but start 
changing for smaller values of k. In particular, we find that for 
a given value of k, £(L) saturates to a finite value at large L. 
This is because, although a large loop occupies a relatively 
large volume at large z, the magnetic energy density very high 
above the disk is very low, and so the contribution of the large- 
z part of a loop to its total energy is relatively small. However, 
we find that the asymptotic value of £(L) at large L begins to 
increase rapidly as k is lowered below roughly K2 — 0.01, in 
agreement with the expectations of § 13.41 

We also calculate the total torque G, transmitted by the 
coronal magnetic field and the total magnetic energy in the 
corona, E tot . We show these quantities as functions of n in 
Figures [13] and [14] Since our code is only first-order accurate 
with respect to the number Ng of grid points in the 9 direction, 
the finite ^-resolution becomes an issue at very small values 
of k, where the distribution function is strongly anisotropic. 
Therefore, to get more accurate values for the torque and the 
magnetic energy corresponding to Ng = oo, we use linear ex- 
trapolation based on the calculations with Ng = 40, 60, and 80, 
separately for each value of k (see Figs. [T3l and [T4l>. 

We find that, as k is decreased, the coronal angular momen- 
tum transfer G{k) increases steadily as k~ 1 for k > 0.1, as is 
expected due to the increased degree of anisotropy of the loop 
distribution function. As k becomes < 0.1, further growth 



of G(k) slows down and becomes G 



-1/3 



The fact that 



G(k) is a power law in this r ange is understandable in light of 
the discussion at the end of § 13.41 combined with the logarith- 
mic dependence of a(n) reported above. We also find that the 
total magnetic energy in the corona, E tot (n), is nearly flat for 
large values of n, but starts to increase roughly as k" 1 / 2 as n 
is lowered below k < 0.05 (see Fig. fT4l>. 

Furthermore, in agreement with the discussion in § 13.41 we 
find that both the magnetic torque and the total magnetic en- 
ergy are substantially larger in the L max = 20 case than in the 
£max = 10 case for small values to k, whereas they are essen- 
tially the same for both values of L max for larger k. 




K = 0.003 

K=0.01 



K = 0.03„ 



0.25 



0.5 



K = 0.1 /' X 
K = 1.0 ^ 
0.75 1 



lo 9 10 z 
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view of convergence with respect to the resolution in the 9 direction for 
k = 0.005 (log| re = —2.3). The plus signs mark the values obtained with 
Ng = 40, 60, and 80 (top to bottom) and the asterisks correspond to the ex- 
trapolation to Ng = oo. 



6. DISCUSSION 

There are several extensions of our model that we plan to 
develop in the near future. They include: (1) characteriz- 
ing the coronal field backreaction on the disk motions; (2) 
including open magnetic flux; (3) taking into account mag- 
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netic twist inside the loops; (4) incorporating a more realis- 
tic prescription for reconnection; (5) investigating the mass 
exchange between the disk and the corona and its effect on 
regulating coronal energy release; (6) studying the interaction 
between the coronal magnetic field and a super-imposed ex- 
ternal large-scale field; (7) assessing the implications of our 
theory for observations, e.g., in terms of time delay between 
hard and soft X-ray emission. We discuss some of these issues 
in more detail in this section. 

6.1. Open Flux Tubes, Outflows, and Net Vertical Flux 

In this paper we have assumed that closed loops are the 
only magnetic structures in the corona. In principle, however, 
one should also consider a population of open flux tubes. If 
considered to be force-free along their entire (infinite) length, 
open loops exert no torques on the disk (in the non-relativistic 
limit) but contribute to the averaged magnetic pressure B 2 (z) 
everywhere and thus prevent it from becoming too small at 
large heights. Open tubes are also essential for investigating 
the role of the corona in launching large-scale outflows; re- 
laxing the assumption of force-free fields allows the inertia 
of the outflow to exert a torque on the disk along open field 
lines felandford & Pavnell 19821) by analogy with the angular- 
momentum loss of the Sun to the solar wind (Weber & Davis 
fT967h . 

In our model, open loops could be introduced through the 
large-scale cut-off for the closed loops. That is, whenever 
a closed loop grows to exceeds a certain maximal footpoint 
separation L max , it could be replaced by a pair of open flux 
tubes of opposite polarities. The large-scale cut-off may be 
physically associated with the disk radius, r. It should deter- 
mine the fraction of magnetic flux that is open at any given 
time. The rules that govern the interaction of open tubes with 
closed ones and with themselves, are straightforward: two 
open tubes of the same sign do not interact; two open tubes 
of opposite sign can annihilate by reconnection, forming one 
single closed loop. An open tube can also reconnect with a 
closed loop, forming again an open tube of the same sign and 
a new closed loop. 

In addition to such pairs of positive and negative open 
flux tubes, there may also be a net vertical magnetic flux 
imposed on the disk by the central object or the ambient 
inter-stellar medium. The radial transport of such externally- 
imposed field across a turbulent accretion disk is an im- 



portant problem with significant consequences for under- 
standing the different spec tral states of accreting black holes 
(Spruit & Uzdensky 2005), for production and collimation of 
disk-driven wi nds and jets, and for s tar-disk magnetic in- 
teraction (e.g. lUzdensky et alJ l2002alfbT) . which is believed 
to regulate the s pin evolution of acc reting neutron stars 
in X- ray puls ars (Ghosh & Lamb 1978), as well as young 
stars dKonigll 1 1 99 U iMatt & Pudritzl 120051) . This problem 
however, is highly non-trivial (e .g. Lubow ^fe SpruiJ 11995b 
Hevvaerts etaTTll996h iLivio et alJ [19991: ISpruit & Uzdenskvl 
2005b . in part because the effective transport of the large-scale 
open flux may be greatly affected by reconnection with the 
small- and intermediate-scale closed c oronal magnetic loops 
(Sprui t & Uzdenskvl 120051 ; lFiskll2005l) . Incorporating a net 
large-scale flux into our statistical model is very straightfor- 
ward and its overall transport should come out automatically. 
We therefore believe that our theoretical model can be a very 
useful tool for addressing this problem. 

6.2. Twisted Loops 

In this paper we have assumed, for simplicity, that coronal 
loops have no longitudinal current, and hence have no inter- 
nal twist. In a more general situation, however, there may be 
force-free field-aligned currents along the loops, generated in 
response to certain disk footpoint motions: specifically, rota- 
tion of field-line footpoints around each other. 

Internal twist would have two major consequences for our 
coronal model. First, the pinch force of the associated longi- 
tudinal current (parallel to the magnetic field because of the 
force-free assumption) will tend to reduce the width of loop 
and its cross section for reconnection with other loops. One 
might consider as a limiting case that the effective thickness 
of the loop is constant along its length. There is indeed obser- 
vational evidence that in the solar corona bright lo ops usually 
have a nearly constant thickness along their length iKlimchukl 
$2(Mi . 

Second, if the loop twist becomes too large, then the en- 
tire loop may become kink-unstable and, as a result, the loop 
makes a transition to a different equilibrium, where the inter- 
nal twist is partially transformed into the global writhe of the 
loop. That is, the loop no longer lies in one plane but rather 
has a twisted, helical shape. Such 5-like loops (so called 
sigmoidal loops) are routinely observed in the solar corona 
dRust & Kumarifl996l) . 

6.3. Incorporating Realistic Reconnection Physics 
(Collisionless Reconnection Condition) 

The physics of reconnection is notoriously complex. How- 
ever, significant progress has been achieved in recent years 
and t he picture that e merges can be summarized as fol- 
lows (lUzden sky 2007b). There are two regimes of reconnec- 
tion: a slow (Sweet-Parker) regime in collisional plasmas and 
a fast regime in collisionless plasmas. The slow reconnec- 
tion regime is just as important as the fast one, since, without 
it, it would be difficult for the system to accumulate signifi- 
cant free magnetic energy before releasing it suddenly via fast 
flare-like events. For practical purposes, the actual rate of fast 
reconnection is not very critical in our problem, as long as it is 
faster than the main dynamical time scale (i.e., the orbital pe- 
riod). More important is the fast-reconnection onset, or trig- 
ger, problem, i.e., the question of the transition from the slow 
to fast reconnection regime occur. The physics of fast colli- 
sionless reconnection is very complicated; it involves either 
two-fluid effects, such as the Hall effect, and/or anomalous 
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resistivity due to current-driven plasma micro-instabilities. 
However, despite this complexity, one can formulate a rough 
criterion for the trans ition from the slow collisional to the fast 
collisionless regim e dCassak et alj|2005t lYamada et alj|2006t 
lUzdens ky 2007J3). We thus plan to utilize this condition to 
formulate a physically-motivated prescription for handling re- 
connection. This prescription can then be used directly in our 
statistical theory or as a sub-grid model in actual MHD simu- 
lations of the corona. In our present model in this paper, the 
reconnection parameter k, which for simplicity we take to be 
constant, effectively subsumes all this complexity. 

The condition for transition to fast collisionless reconnec- 
tion involves several physical parameters of the system, in- 
cluding the ambient plasma density. The dependence of the 
reconnection regime on the density is critical, since it es- 
tablishes an important feedback that the dynamically sub- 
dominant coronal gas exerts on the coronal magnetic field 
dUzdenskvl 2007a b). In turn, the plasma density in the 
corona is determined by the disk-corona mass exchange pro- 
cesses, such as evaporation in response to coronal heating, 
precipitation due to gradual radiative cooling, and magneto- 
centrifugally and radiatively driven winds. Coupled together, 
disk-corona mass exchange and the transition to fast colli- 
sionless reconnection ensure that the corona is maintained 
near the state of marginal collisionality and regulate the over- 
all level of coronal activity and its intermittency, as well as the 
vertical distribution of magnetic energy density and of mag- 
netic dissipation. 

These ideas, especially the concept of marginal collisional- 
ity, have recently been successfully applied to the solar corona 
to expla in the self-regula ting nature of the coronal heating 
process dUzdenskv|20 07 a b) and to the coronae of other main- 
sequence stars dCassak et al. 2007). They also have proved 
very useful for providing a natural explanation for the ob- 
served optical depth in the coronae of accreting black holes 
(Goodman & Uzdensky 2008, in preparation). 

Finally, especially for studie s of the ADC's interaction with 
a larg e-scale disk wind (e.g. Brandenburg & von Rekowskil 
2007), it is important to understand the transition from the 
force-free regime to the wind regime in which the plasma iner - 
tia becomes dynamically important dUzdensky et alJl2002bh . 
This transition happens near the Alfven critical surface of the 
outflow and has a strong effect on magnetic reconnection. In 
particular, it is expected that, beyond the Alfven surface, open 
magnetic fie ld fines will n ot be able to close back via recon- 
nection (e.g. lUzdens kv 2004). 

6.4. Flux Emergence 

Our model, by construction, is not complete — it needs to 
be connected to what happens in the disk. In particular, it 
needs as input a statistical description of the population of 
small loops, or of their rate of emergence from the disk. This 
information determines the overall normalization of the loop 
distribution function and hence the net rates of angular mo- 
mentum transport and dissipation in the corona. In the mod- 
els computed here, we have simply assumed a fixed, isotropic 
distribution of small loops. 

The ultimate source o f coronal activity i s the MHD tur- 
bulen c e in the disk (e.g. iGaleevetal] Il979t iTout & P ringle 
Il992t iMiller & Stonel |2000|). Therefore, in order to esti- 
mate the rate and form of magnetic flux emergence into 
the corona, one first needs to understand the properties of 
MHD turbulence in a stratified disk, with a particular em- 
phasis on the production, evolution, and buoyant rise of mag- 



netic flux tubes (e.g. lSchramkowski & Achtergergl 19931) . The 
best prospect for developing this understanding is through ap- 
propriate statistical analyses of MRI turbulence in stratified 
shearing boxes. 

However, the notion of a flux tube in most theoretical 
discussions (including the present paper) is not sufficiently 
precise to be applied directly to simulations. What well- 
defined, measurable, and statistically meaningful quantities 
correspond to flux tubes or to their rates of emergence? Can- 
didates exist, but it is not clear which is best. Field lines can 
be found as integral curves of the field and their motions deter- 
mined, but by what prescription should they be grouped into 
tubes? Alternatively one might work with Fourier decompo- 
sitions of the vertical Poynting flux on horizontal planes, or 
with more general n-point correlation functions of the field. 

7. SUMMARY AND CONCLUSIONS 

In this paper we construct a general theoretical frame- 
work for understanding the structure of a strongly-magnetized 
corona above a turbulent accretion disk. This study is moti- 
vated by the need to provide a more solid physical foundation 
for ADC spectral modeling efforts. It should also act as a con- 
necting bridge between numerical MHD simulations of MRI- 
turbulent disks and semi-empirical coronal models, and stim- 
ulate further theoretical studies of accretion disks coupled to 
their coronae. We also hope that some of the theoretical tools 
and ideas developed in this paper will prove useful in solar 
physics. 

One of the major goals of our study is to develop a statis- 
tical language appropriate for describing the chaotic coronal 
magnetic field. Here, we are interested in spatial scales larger 
than the disk thickness but smaller than its radius and in tem- 
poral scales longer than the orbital period but shorter than the 
overall accretion time. Our approach builds upon the previous 
work by Tout & Pringle (1996), but uses much more realistic 
physics in several key aspects and also goes further in analyz- 
ing and interpreting the results. 

To construct the statistical theory, we represent the corona 
by an ensemble of elementary magnetic structures, namely, 
loops connecting two spots on the disk surface (Fig. [TJ. Each 
loop is characterized by several primary parameters (e.g., the 
distance between the footpoints and the orientation). The 
main object in this study is the distribution function F of loops 
in this parameter space. One of our main goals is to formulate, 
and then solve, the loop kinetic equation (LKE) for this distri- 
bution function, similar to the Boltzmann kinetic equation in 
the statistical theory of gases. 

To do this, we first analyze the key physical processes that 
govern the evolution of individual coronal loops. First, there 
are several processes that affect the loops individually, such 
as: (1) emergence of small loops into corona; (2) random 
footpoint motions due to the disk turbulence; (3) Keplerian 
shear, stretching loops azimuthally and thereby also making 
them grow in height. On average, these processes pump en- 
ergy from the disk into the corona, creating a stressed non- 
potential force-free field. In addition, loops may interact with 
each othe r via episodic reconn ection between two individual 
loops (c.f. Tout & Pringle 1996), forming two new loops (see 
Fig. [5]). In our theory reconnection is represented as a bi- 
nary collision, analogous to binary collisions between atoms 
in a gas. Correspondingly, we describe loop-loop reconnec- 
tion by a nonlinear integral operator, similar to Boltzmann's 
collision operator. In contrast to processes (l)-(3), magnetic 
reconnection relaxes the accumulated magnetic stresses and 
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dissipates the accumulated free magnetic energy. Overall, a 
magnetically-active ADC can be described as a Boiling Mag- 
netic Foam. 

Based on these processes we are able to construct the loop 
kinetic equation. In this equation we characterize the over- 
all rate of reconnection events relative to the Keplerian shear 
rate by a dimensionless parameter n. In order to investigate 
the role of magnetic reconnection in the corona, we solve the 
loop kinetic equation numerically for several different values 
of k. We obtain a statistical steady state for each value of k 
and find that the steady-state loop distribution function is gen- 
erally well represented by a orientation-dependent power-law, 
F(L,8) ~ L" Qk(8) . When Keplerian shear is absent, the distri- 
bution function is isotropic, = const. As the rate of 
shear relative to reconnection increases (i.e., n decreases), the 
distribution becomes more and more anisotropic, with a pre- 
dominance of toroidal orientation. At the same time, a typi- 
cal loop grows to a larger size by a stronger shear before its 
growth is disrupted by reconnection. Thus, the orientation- 
averaged distribution function becomes shallower as k is de- 
creased. 

Once the distribution function is known, we use it to calcu- 
late several important integral quantities related to the ener- 
getics of the corona. First, we use a self-consistent mean-field 
approach to compute the magnetic energy density as a func- 
tion of height, B 2 (z)/87r. This quantity represents the collec- 
tive magnetic pressure of all the neighboring loops that con- 
fine any given loop and thus represents another (in addition 
to reconnection) way in which loops interact with each other 
in our theory. Although it doesn't enter explicitly into the 
loop kinetic equation, B(z) is very important in our model. In 
particular, it controls the thickness of loops as a function of 
height, which affects in turn the cross-section for reconnec- 
tion. In addition, by requiring that the vertical gradient of 
the magnetic pressure B 2 (z)/8n be balanced by the magnetic 
tension within each loop, we self-consistently calculate the 
equilibrium shape and vertical extent, Z top (L), of the loops. 
This, in turn, enables us to calculate some important quanti- 
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ties such as the energy associated with a given loop, the force 
exerted on its footpoints, etc. We then use these quantities to 
assess various issues of coronal energetics, including the over- 
all magnetic energy stored in the corona, statistical distribu- 
tion of coronal energy release events (flares), and the overall 
angular momentum transferred by the coronal magnetic field. 
As a result of our parametric study with respect to the recon- 
nection parameter k, we find that if k is decreased (i.e., re- 
connection in the corona is inhibited) beyond a certain value, 
the slope of the loop distribution function becomes so shallow 
(namely, shallower than L~ 3 / 2 ) that the contribution of large 
loops to both the magnetic energy and torque starts to domi- 
nate, leading to a significant enhancement in these quantities. 
In our specific model, the critical value of n is found to be 

Ky 2 - 0.002. 

These results demonstrate that the energetic dominance 
of coronae is inextricably linked to reconnection processes. 
They thus motivate further efforts to develop more realis- 
tic physical description of reconnection. To reiterate an im- 
portant point made in the Introduction, the tenuous corona 
above an accretion disk is likely to be marginally collision- 
less (Goodman & Uzdensky 2008, in preparation), unlike the 
dense plasma inside the disk itself. This means that the corona 
cannot be described by traditional MHD simulations with 
constant resistivity because of their inability to control or re- 
solve magnetic reconnection, which, as we have shown in this 
paper, may influence the coronal magnetic energy and angu- 
lar momentum transfer. Therefore, some kind of a physically- 
motivated subgrid prescription for reconnection is needed. 
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APPENDIX 

APPENDIX A: CALCULATION OF B IN TERMS OF F(L) 

In this Appendix we calculate the mean magnetic field B(L) in terms of the loop distribution function F(L) corresponding to a 
self-consistent atmosphere (see § I3.21 >. 

Consider a horizontal slab of thickness dz at height z (see Fig. |A15t . Correspondingly, there is a minimal length L of loops 
that reach above this height. Consider now an arbitrary slender loop A of length L' > L, with a cross-sectional area a(z) and 
the angle between the loop's direction and the vertical equal to a(z,L'). Then the volume that the two segments (ascending and 
descending) of this loop occupy inside the slab under consideration can be written as 

dV = 2a(z) (Al) 

cos a(z,L') 

The volume occupied by all such loops per unit horizontal area of the slab is equal to 

oo 

dV = 2 { dU a{z)F(L') . (A2) 

J cosa(z,L') 

L 

Since in our model each loop carries the same amount A^ of magnetic flux, the loop cross-sectional area a(z) is the same for all 
loops at a given height and is simply equal to A^>/B(z). Thus, the above volume is 

oo 

2A* f - , dL' 

dV=^—dz F(L') — — . (A3) 

B(z) J cosa(z,L') 

L 

But, on the other hand, these loops occupy all the volume within the slab, and so this volume per unit horizontal area ought to 
be equal to just dz- This gives us the equation that determines B as a function of L: 

oo 

B(z) = B b(z) = 2A* f F(L') = . (A4) 

J y/l- sin 2 a(z,L>) 
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Now, using equation ( [Pil l, we can replace sin a(z,Z/) = b top (L') / b(z) = b top (L') / b tov (L) and hence obtain the following integral 
equation for the function b(L'): 

oo _ 

F(L')dU B 



L V "topVW "top 

which can be rewritten as 



bUL)-bUL>) 2Av& 



: const , (A5) 



* U(b')db' B 



*Jb 2 -b' 2 2A* 



: const (A6) 



where U(b') = -F[L(b')]dL/db' . By substitutions t = b' 2 , s = b 2 , and V(f) = U(b')/2b', this equation can be transformed into the 
Abel integral equation, that can be immediately solved, yielding the following final result: 

V(s)= 2Jk^ => ^=-F [Lmd L/ db =^- o . (A7) 



db = - F(L) dL => b(L) = I F(L') dL' , ( A8) 



and hence 

dh = — 

Bo Bq 

L 

which is in agreement with our general expectation above. 

APPENDIX B: PROOF THAT 8{L) = 2E MACN 
In this Appendix we prove the conjecture that £{L) = 2E m . d „ a (L). The proof goes as follows. 

First, according to equation A3 It , the force on the loop's footpoint can be expressed in terms of the magnetic field at the top of 
the loop as 

ML)= ^ = 47T • (B1) 

Using expression (l34t and taking into account that a loop has two legs, we have 

a* r a* r - dz a^ r b 2 (z) 

£ magn (L) = 2x — / B(l)dl=— / B(z) = — / W dz. (B2) 

87r J 4tt J cosa 4ir J r2_ r 2 

left leg z=0 z=0 V P 

On the other hand, substituting (IB lb into equation (l30l l. we get 

L 

£(L)=^^ I b top (L')dL' . (B3) 

Integrating this by parts yields 



L 1 

I\jtI u/t\t i / TtU\jU 



£(L) = j b(L')dL' = b(L)L+ J L(b')db' . (B4) 

A(L) 

According to (fTTT i: 

L/ /■ W')M>"] ,,,,,, 



&'=&(£) 

Interchanging the order of integration, we get 

l l 



L(b')db'=-2 / V /; ' Ah" Ah' (B5) 

' ' Vb" 2 -b' 2 



L(b')db' = -2 I db" d -^-l ( [ b ' Ah' 
J db" \J ^b" 2 -b' 2 

b"=b b'=b 
1 

= _ 2 / J^i—tf^Lldb". (B6) 

./ db" 



b"=b 

On the other hand, according to ( fTTI i. 



b(L)L = -2 / W . (B7) 

7 Vb" 2 -b 2 

b"=b(L) 
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Combining these results we get 
End of Proof. 



1 

f b" 2 dz" 

£(L) = -2 J 7 ===2E m > sn . (B8) 

b"=b(L) 



